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1. Introduction 

These lectures provide an introduction to the theory of disordered interacting electron 
systems. As for the case of superconductivity, the understanding of the behavior of 
transport and thermodynamical properties of metals and semiconductors has required 
the invention of new and fascinating concepts. 

In these lectures our focus is mainly on the effects of disorder, and its interplay with 
electron-electron interaction. The resulting theory, although still cannot answer some im- 
portant questions, is simple and elegant. It describes the combined effects of interaction 
and disorder in terms of a renormalized Fermi liquid, whose Landau parameters become 
scale dependent and provide, together with the conductance, the couplings flowing under 
the action of the renormalization group. However, this final simple description, which 
has required several decades of intensive work from many people, is built on several con- 
ceptual steps. It is our aim to lead the reader through the development of these various 
steps. Our hope is that the reading of these lectures should allow people not expert in 
the field to access the original literature. 

There are already several review articles which give an account of the problem from 
different view points and at different stages of the historical development [21 ED EH 03 El 
0|H]. However, the most recent and complete are still quite hard to read for unprepared 
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readers. For unprepared readers we mean those people that are not familiar with the 
complex technical jargon that the field experts have developed over the years. 

We will concentrate on those aspects that we believe are fundamental for the problem 
of the metal-insulator transition due to disorder and interaction. This will force us to 
ignore a number of extensions and developments of the theory. These latter, however, 
may be found in the existing reviews. 

We have also chosen to present the theory in the simple language of standard many- 
body perturbation theory. The field-theoretic approach based on the derivation of an 
effective non-linear <7-model[0] is certainly more elegant and powerful, but requires quite 
some effort to appreciate the beauty of it. We invite the reader to approach it after 
reading these lectures. These lecture notes are self-contained. A basic knowledge of 
many-body theory and diagrammatic technique is the only prerequisite. 

After these warnings, we outline the contents of these lectures. In the next section we 
set the stage for the microscopic theory by introducing the reader to the metal-insulator 
transition in disordered systems and to phenomenological scaling. The necessary back- 
ground for the microscopic theory is given in the following section. The fourth section 
deals with the non-interacting problem. A number of key physical and technical in- 
gredients are introduced in a pedagogical way. Also, the experimental urgency to take 
into account interaction effects is presented. The fifth section goes to the heart of the 
problem by building the renormalized disordered Fermi liquid. Gauge invariance and 
Ward identities are the shining lighthouses which help us to navigate through the messy 
waves of the perturbation theory. Land is finally reached in the sixth section, where we 
discuss the renormalization group equations and look back to our journey and compare 
the theoretical understanding with the available experiments. 



2. Setting the stage for the metal-insulator transition 

In this section, we begin by recalling the textbook theory of electrical transport in 
metals. Then we move to a description of the actual physical systems where the phe- 
nomena, which we describe theoretically, are observed. We conclude the section with the 
scaling theory of the metal-insulator transition due to disorder. 

2T. The semiclassical approach of Drude-Boltzmann. - The conventional theory of 
electrical transport is due to Drude. In its original formulation, Drude suggested that 
electrons, under the action of an externally applied electric field, are accelerated according 
to Newton's equation of motion until they collide with the ions after a time r. The 
distance between successive collisions determines the mean free path I. Due to this 
sequence of independent scattering events the electrical conductivity is given by 



(2.1) 



c = 



e 2 noT 



where no is the density of electrons and e, m are the charge and electron mass. 
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Fig. 1. - A pictorial representation of the semi-classical theory of transport. 



After the birth of quantum mechanics, Sommerfeld reformulated Drude's theory to 
accomodate the Fermi statistics of electrons, providing the correct relation between t and 
I via the Fermi velocity Vp. More importantly, with the work of Bloch, it was realized 
that the relaxation of electron momentum and the finite value of the conductivity is due 
to imperfections of the ion lattice, i.e. to disorder. Drude's law l|2.1l) may be obtained 
by a semi-classical approach based on the Boltzmann equation for the evolution of the 
electron distribution function in the presence of external fields. A pictorial description 
of Drude's model of electrical conduction is shown in figd 

As a consequence of the collisions, the electrons undergo a classical random walk of 
step I and a diffusive motion, with the diffusion coefficient D related to the conductivity 
by Einstein's relation 



(2.2) 



OfJ, 



where dno/d/j, is the compressibility. In the case of the Fermi gas, dno/d/j, = 2Nq is 
simply related to the density of states per unit volume per spin 



(2.3) 



Qrf d-2 j /0 d ~ 2 



where Qd is the solid angle in d dimensions and Ep the Fermi energy. From the Drude 
formula (|2.1|l and Einstein's relation (|2.2() . with uq = 2NolEp/d, one gets the diffusion 
coefficient!? = (2Epr)/ (d m) = vpr/d. Within the independent electron approximation, 
only one diffusion constant D controls the charge, spin and heat transport, leading to 
relations similar to ea. (|2.2|) for the charge. In particular, the thermal conductivity, ke, 



(2.4) ke = Cv.qD 

where Cv.o — (2tt 2 /3)kgNoT is the specific heat for the electron gas. 
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In concluding this subsection, we introduce the conductance related to the conduc- 
tivity by geometrical factors 

(2.5) G=aJ T = aoLd ~ 2 

where S and C are the cross section and length of the conductor to which we assign the 
typical size L. By using the explicit expression of the density of states, one may rewrite 
G as 



(2.6) G 



2e 2 Oi 2E f t fp F L\ d 2 _ 2e 2 ft d p F lfp F L\ d2 



which shows that, in the natural conductance units (Go = 2e 2 /h — 12.9kil~ 1 ), the 
value of G is controlled by the dimensionless parameters p F l/h — 2irl/\ F and L/Xp. 
In two dimensions, in particular, conductivity and conductance have the same physical 
dimensions and the ratio between the Fermi wavelength and the mean free path is the only 
parameter that controls the value of the conductivity. In the semi-classical limit, X F <C I, 
the Drude formula predicts a high conductivity. The rate of collisions r _1 is proportional 
to the impurity concentration. By increasing the disorder in the semiclassical approach, 
one has that <7o diminishes, but remains finite. Ioffe and Regelj9| stated the criterion that 
in the metallic phase the mean free path / cannot be smaller than the average interelectron 
distance proportional to h~/p F , i.e., p F l/Ti > 1. Mott|ll)| applied this criterion to the 
Drude conductivity arguing that for d > 2 there is a minimum metallic conductivity, 
CTomin when / w h/p F , 



(2.7) <7 0n 



2e 2 Q d f PF L\ d - 2 
h (2 7 r) d - 1 d V h J 



As a result, there should be a discontinuity of the conductivity (which is universal in 
d = 2) at the transition from the metal to the insulator. However, when I « h/p F we 
are deeply in the quantum limit and the Ioffe-Regel criterion cannot be naively applied 
to the semiclassical Drude formula. Indeed one expects that corrections beyond the 
semi-classical approximations will strongly modify ea. l|2.1|l opening the way to a new 
perspective in the metal-insulator transition. Most of these lectures concern precisely 
this type of corrections. 

2'2. The metal-insulator transition. - There are, of course, finite-temperature cor- 
rections to Drude's formula. In general, temperature-dependent corrections arise from 
inelastic scattering of electrons beween them and with the phonons, resulting in the 
characteristic w T 2 and w T 3 behavior of the conductivity. However, typical disordered 
systems show, at low temperature, strong anomalies. In metallic films, for instance, there 
are temperature dependent logarithmic corrections 



(2.8) 



<t(T) = do + m In T, 
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Fig. 2. - Energy diagram of an n-type semiconductor containing donors and acceptors. The 
horizontal lines represent centres, the circles electrons in them. 



where m is positive. Besides metals, experimental realizations of disordered systems 
are obtained in doped semiconductors like Si : P, Ge : Sb and amorphous alloys as 
Nb : Si, Al : Ge, Au : Ge. In a doped semiconductor, there are two types of conduction 
mechanisms. The first is due to the thermal activated carriers and dominates at high 
temperature. To understand the second mechanism, let us consider, for instance, an 
n-type semicondutor, as shown in fig[2 An electron sitting at its donor atom location, 
has a wave function exponentially localized around the impurity (the energy of such an 
electron is indicated by a short horizontal line in the figure) . Due to the small, but finite 
overlap of wave functions centered at different impurity locations, the donor electrons 
can move around by tunneling from one impurity to another. This gives rise to what is 
called impurity conduction. A doped semiconductor is compensated when, besides the 
majority donor atoms, it contains also some minority acceptor atoms. In this way, some 
of the donor electrons are captured by the acceptor levels, by allowing the tunneling of a 
donor electron from an occupied level to an unoccupied one. By increasing the impurity 
concentration, the overlap of the wave functions sitting at different impurity sites becomes 
larger. One point to notice is that by increasing the impurity concentration, there are 
two competing effects. On the one hand, disorder increases due to the larger number 
of scattering centers. On the other, at a high enough concentration of impurities, the 
overlap is such that the impurity levels form a band, which behaves as an intrinsically 
disordered degenerate electron gas and yields a metallic conductivity. Hence disorder 
effects are stronger at lower concentration. The transition to metallic behavior of the 
impurity conduction occurs at a critical impurity concentration, n c . Si : P, where P 
donors sit substitutional^ and randomly in a dislocation-free Si lattice, is an ideal system 
to study the effect of disorder on transport properties. For instance, at enough impurity 
concentration to be in the metallic state, one measures 



(2.9) 



a(T) = a + mT\ 
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where the coefficient m can be both positive and negative and n = 1/2 J2]. By decreas- 
ing the P concentration below a critical value, the system undergoes a metal-insulator 
transition at T — in the sense that 



with the critical exponent /i = 1/2 ^|^] The value of n for uncompensated Si : P is 
still under debate and depends strongly on the identification of the critical region in the 
experimental data[H[TSl[TBl[nHTl]. Compensated samples [T5| and the alloys [2011^111^] 

have fx = 1. 

Besides the transport properties, anomalies are also seen in the tunneling density of 
states for Au : Ge|2l]. NbSifXft, in specific heat 12111311201113, and in spin susceptibility 
|281 1291 15U1 1311 132| in SiP. As we will discuss in a more detailed way in the next sections, 
taking into account corrections both in transport and thermodynamic properties will be 
crucial in developing an effective Fermi-liquid theory for these systems. 

In more recent years, the discovery of a metal-insulator transition in the two-dimensional 
electron gas 33 has stimulated a renovated effort to understand the interplay of disorder 
and interaction effects This phenomenon has been first observed in Sz-MOSFET 
devices and later also in other two-dimensional electron gas realizations as in semicon- 
ductor hetero-structures. In Sz-MOSFETs devices, the two-dimensional electron gas is 
formed at the interface between the bulk silicon and an insulating layer of silicon oxide, 
as shown in figOl By applying a positive bias on the metallic gate deposited on the 
insulating silicon oxide layer one forces the electrons to move in the Si02-Si interface, 
in an almost two-dimensional environment. As compared with the doped semiconductor 
systems, these systems have the advantage that the density of the electron gas is almost 
continuously controlled by the degree of the band bending at the interface, i.e., by the 
applied bias. This allows a very fine scanning of the properties of the sample as function 
of density. Furthermore, the disorder is mainly due to scattering centers in the insulating 
layer, so that in principle one varies the density at fixed amount of disorder. By varying 
the electron density, n, one can change the effect of the interaction since Ep oc n, and 
in 2D the Coulomb electron-electron interaction Eq oc n 1 / 2 . The ratio r s between the 
Coulomb interaction evaluated at the average interparticle distance and the Fermi energy 
is given by 



where e is the dielectric constant. By using Ep = nj (2 No) and r av — l/\/n and recalling 



(2.10) 



cr ~ (n - n c ) M , 



(2.11) 



E c _ e 2 /(sr av ) 



Ekin Ep 



( 1 ) At present there is not yet a general consensus on whether we have a real zero-temperature 
transition or rather a crossover effect. 
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Fig. 3. - Scheme of a SiMOSFET device. The metallic gate (Al) is positevely biased so that it 
attracts the electrons, which on the other hand cannot enter the insulating layer of SiC>2- The 
electrons are then confined at the interface Si — SiC>2 and form a two-dimensional electron gas. 
In the proximity of the interface the slope of the energy bands determines the effective tickness 
of the two-dimensional electron gas. 



the expression of the Bohr radius as = eh 2 /(me 2 ), one obtains 

(2.12) r s = X —=. 

TreaBVn 

At the present, the origin of the metal- insulator transition in two-dimensional systems is 
still unclear and represents a very hot issue of debate in the literature. (A recent review 
may be found in refs.|j21 [35]). For this reason, we prefer not to enter now in a detailed 
discussion of the experimental features of this phenomenon, which we will point out later 
on when relevant results of the theory will require it. 

2'3. The Anderson transition and quantum interference. - All of this suggests that 
disorder cannot be treated only within a semi-classical approach. In 1958, Anderson 
invented the field of localization, proposing that under certain circumstances diffusion 
may be completely suppressed 36 . He proposed a lattice model where the site energies 
are randomly distributed. When disorder is absent, a small hopping amplitude is enough 
to delocalize the electron states and form Bloch waves. However, by increasing the 
disorder, the hopping processes may only spread an initially localized state over a finite 
distance, which defines the localization length £o- Since in the process of the impurity 
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band formation, localized states are more likely to form in the band tails, Mott argued [3"?) 
that there must exist a critical energy E c , called the mobility edge, which separates 
localized from extendend states. When the Fermi energy Ep is below the mobility edge, 
the system is an insulator. When Ep passes through E c , the system becomes metallic. 
From this point of view, for a given model and given amount of disorder, the problem is 
to compute E c . For non-interacting systems, this can be tackled numerically, by exactly 
solving the Schrodinger equation in a disordered lattice. A review of the status of the 
numerical simulations may be found in rcf. 7 . 

Even though these concepts played an important role in shaping our modern view 
of the metal-insulator transition, a great impulse to the development of the field came, 
however, by the discovery of the phenomenon of weak localization. 

As remarked at the end of the previous subsection, the standard theory of transport 
is based on a semi-classical approach, where in evaluating the probability for electron 
diffusion one neglects the interference between the amplitudes corresponding to different 
trajectories and essentially treats a classical random walk with step I and diffusion con- 
stant D. This is indeed justified in many cases, where the semi-classical theory works. 
In fact, in a disorderd system the phase difference for any two different trajectories will 
vary randomly. The situation changes, however, for self- intersecting trajectories, which 
come from closed loops. In this case, trajectories naturally come in pairs, depending 
on whether one goes around the loop clockwise or counter-clockwise. One expects that 
interference between these pairs of trajectories modifies the semi-classical result. There 
is a simple argument to estimate the probability of having a self-intersecting trajectory. 
On the one hand, one has that the electron motion is described by a classical diffusion 
process such that the average distance after a time t is 



On the other hand, the quantum nature of the electron may be thought of in terms of a 
tube of size Xp generated by the electron motion. In a time di, the volume spanned by 
the tube increases by 



Let us consider the ratio between the increase of the tube volume in time dt and the total 
volume generated by the diffusion process. The total probability for self-intersection may 
be estimated by integrating this ratio over time 



where the lower limit r is the time above which the diffusive regime, after a few collisions, 
starts to set in. The upper limit, r^, is the time until which phase coherence of the wave 
function persists. In general, inelastic processes at finite temperature make 70 a decreas- 
ing function of temperature. In two dimensions the probability grows logarithmically as 



(2.13) 



Dt. 



(2.14) 



dVtube = X F 1 v F dt. 



(2.15) 
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temperature decreases. At zero temperature, when — > oo, the upper limit is provided 
by the system size via the diffusion relation L 2 /D = tl- In this way, the probability of 
self-intersection acquires a scale dependence 



(2.16) P 



v F X° 



D d/2 T (d-2)/2 d _ 2 



Go 



a l d - 2 d-2 



By assuming that the conductivity corrections are proportional to this probability, one 
obtains 



(2.17) 



5a 1 



where e = d — 2 and go = G(1)/Gq is the conductance at the scale I in units of Go- 
Equation H2.17fl is valid at T — and the inverse scattering time or the inverse mean free 
path play the role of the ultraviolet cutoff, whereas 1 is the infrared cutoff. At d = 2, 
the conductivity correction is log-singular. 

At finite temperature, when < tl, the infrared cutoff becomes temperature depen- 
dent and in 2D the correction becomes logarithmic in temperature. This opens the way 
to the scaling theory discussed in the next subsection. 

As a final remark to this subsection, we point out that the weak-localization phe- 
nomenon is sensitive to any perturbation that breaks the time reversal invariance. This 
is clear from the above argument of the interference between time reversed trajectories. 
For instance, in the presence of a magnetic field, the two trajectories acquire a phase 
difference 4>\ — 4>2 — (2e/ hc)^ b , proportional to the magnetic flux <&b threading the sur- 
face delimited by the closed loop. Since at finite temperature the logarithmic singularity 
is cut off at time T0, in order to cut off the singularity typical magnetic fields must be 
of the order of a flux quantum over a region whose size is of the order of the dephasing 
length = y/Dr^. This gives the condition B > (h / ec) / L^. It is also clear that further 
dephasing mechanisms, as for instance, spin-flip scattering with typical time t s become 
important when t s < t^. 

2'4. The scaling theory of the metal-insulator transition. - The starting point is the 
argument of Thouless concerning the evolution of the wave function as the system size 
is increased [HH1 EH] • To fix the ideas, let us imagine that the system of system size 2L 
is made up by combining blocks of size L, as shown in fig0] Suppose we know the 
eigenstates for a block of size L. The level spacing for these is AE. We ask how are 
the states when we combine blocks together. Let SE be the energy brought about by 
the perturbation of joining the blocks together. Clearly, we expect that for 8E <C AE 
the new eigenstates for a system of size 2L will differ very little from those at scale 
L. The energy 5E measures the sensitivity to a change in the boundary conditions. 
In a diffusive sytem this may be related to the time necessary to reach the boundary, 
5E = hD/L 2 . The level spacing, on the other hand, is related to the inverse density of 
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Fig. 4. - Block of size 2L obtained by blocks of size L. 



states, AE = l/(NoL d ). By using Einstein's relation, the ratio of the two energies gives 



(2.18) 



5E 
~AE 



NnL d 



hD 



i <y L d - 2 

2ir 2e 2 /h 



g(L) 

2tt 



where g{L) is the conductance at the scale L in units of 2e 2 /h. If g <C 1, the new 
eigenstates are not modified much by the assembling of the blocks. On the other hand, 
when j > 1, the new states are delocalized on all the blocks. The scaling theory [ID] 
assumes the g(L) is the only parameter that controls the evolution of the eigenstates 
when we rescale the system size. Mathematically this is expressed by requiring that the 
conductance of a block of size L' = bL is expressed in terms of the conductance of a block 
of size L by a function of L' /L and g(L) only, i.e., g(L') — f(L'/L,g(L)). Its logarithmic 
derivative for V = L, which defines the /3-function of the corresponding renormalization 
group equations, depends on the scale L only through g(L) itself 



(2.19) 



d In g(L) 
dlni 



= P{9{L)). 



The vanishing of the /3-function controls the scale-invariant limit, i.e., provides the fixed 
point of the trasformation g* . In the case of one-parameter equation the fixed g* point 
coincides with the critical value g c . Linearization of the transformation, starting from 
the fixed point, provides the scaling behavior of the physical quantities. The /3-function 
is relatively well known in the two limits of a good metal, where Ohm's law is valid and a 
is a constant, and in the strongly localized insulating regime, where the scale-dependent 
conductance falls off exponentially over a localization length £o as 



(2.20) 



g{L) = g a e 



-LI to 
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Fig. 5. - Schematic /3-function. 



One then immediately gets 

(2.21) /3(g) =d- 2, 5 »1 

(2.22) (3(g) = In-?-, 5 «1 

90 

where go is the conductance at some initial microscopic scale I. Under reasonable assump- 
tions, the (3- function has the qualitative behavior shown in fig|SJ A positive (negative) 
value for (3 means that upon increasing the system size, g increases (decreases) corre- 
sponding to a metallic (insulating) behavior. 

A zero g c such that f3(g c ) = 0, signals an unstable fixed point for the flow of g. This 
represents a metal-insulator transition. One consequence of eas. l|2.21|) . H2.22|) is that, for 
d < 2, the system is always an insulator at zero temperature and all states are localized. 
Close to a critical point at d > 2 we may linearize the /3-function to get 



(2.23) 



from which 



d(g 



dlnL V d -9 ) g =g c 



(2.24) g(L) -g c = (g - g c ) (j \ " 



12 



Carlo Di Castro Roberto Raimondi 



where go is the conductance at scale I and 

(2.25) x g = 5 c/5'(5c)- 

By defining a correlation length, which coincides with the localization length in the 
insulating side, it diverges at criticality as 

(2.26) e ~ (fld ~ 5c)-". 

By assuming that £ is the only relevant length, it scales as £' = £/(L/l) and one deduces 
^ = 1/xg. Furthermore, from the critical behavior of the conductivity 

(2-27) ,„to_ flo) /.„ |[g_>_*_ 

one derives the scaling law|41j 

(2.28) (j, = (d - 2)v = ev. 

In the metallic regime, where g is large, we can assume that the /3-function can be 
expandend in a power series in l/q |42ll4"0] 

(2.29) 0(g) = d-2---\ + ... . 

9 9 

Above two dimensions, if a > 0, one has a fixed point g c — a/e and v = 1/e leading to 
[i=l. On the other hand, if a = 0, the fixed point is determined by the second order 
term, g 2 c = b/e, which implies v = l/(2e) and /i = 1/2. Finally, when a < 0, there is no 
fixed point at this order. 

At d = 2 the scaling equation reduces to 

do a 

^ dhrL = -g 

For a > the system scales to an insulator, whereas for a < it scales to a perfect 
conductor. This phenomenological theory does not allow for a metallic phase in d = 2. 

3. The microscopic approach 

In this section we introduce a few general tools that will be used in building a mi- 
croscopic theory. First we discuss how physical observables may be evaluated in terms 
of response functions. Secondly, we show how conservation laws impose constraints, 
e.g. Ward Identities, on these correlation functions, which are useful when performing 
perturbative expansions. We mainly follow ref. 43 . 
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3'1. Linear response theory, Kubo formula and all that. - The coupling with an ex- 
ternal electromagnetic field is given by the Hamiltonian( 2 ) 

(3.1) H emi = - c J drA"(r)J M (r) 

where the Greek index runs over time (/i = 0) and space indices (fj, = 1, d). The latter 
will be later on indicated by Latin letters . As it is standard, lower indices have space 
components with a minus sign, e.g., J M = (cp, — J). With r we indicate the position 
vector in any dimension d. The external scalar, </>(r), and vector potential, A(r), are 
coupled with the charge and current density, defined by 

(3.2a) p(r) = eV4(r)^(r) 

(3.26) J(r) = [^t( r )VW(r) - (V^(r))^(r)] - ^-A(r)^(r)^(r) 

(3.2c) =j(r)-— A(r)^(r)^ CT (r). 

mc 

In the following, for the sake of simplicity we shall set h, c, and ks equal to one. In 
eas. (|3.2ll . t/v( r ) (V4-( r )) is the annihilation (creation) Fermion field operator. Our goal is 
to study the system response to an external electromagnetic field within linear response. 
The second term in ea. <|3.2c[) . the diamagnctic contribution, being already linear in the 
field, may be evaluated as 

(3.3) ] dia = --n A 

TO 

where Hq is the equilibrium (number) density. By the compact notation x = (i, r), the 
linear response is given by 

(3.4) J"(a:)= [ dx' K^(x,x')A v (x') 



where the response kernel K^ v (x,x') is the four-current correlation function, which in- 
cludes both the density-density and current-current correlation functions, 

(3.5) K^(x, x') = R^(x, x') + —n S {i) (x - x')S^(l - S M ) 

TO 

and 

(3.6) R» v {x,x')=-i6(t-t') < \j»(x),f(x')} >, 

( 2 ) We adopt the relativistic notation: upper and lower indices indicate contravariant and 
covariant vectors, respectively, i.e. A M = (0, A) and — ((f),— A). 
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with the average taken over the appropriate statistical cnscmblc( 3 ) . If the unperturbed 
system is traslationally invariant and has a time-independent Hamiltonian, we can use 
Fourier transforms with respect to both r — r' and t — t' ', 

/DO 1 
-co 2tt q 

The sum over the momenta is left unspecified for the time being. It depends on the 
choice of boundary conditions. In the limit of an infinite system, the sum gets replaced 
by an integral over all space in the standard way. In Fourier space, ea. (l3.4|) becomes 
local 

(3.8) J*(q) = K^(q)A u (q), 

where q — (w,q). Physical observables are now readily obtained. For instance, the DC 
electrical conductivity, by making the choice of a time-dependent vector gauge, E = 
— d t A(t), reads 

(3.9) fftf = -lim 

3'2. Conservation laws and gauge invariance. - Charge conservation is expressed by 
the continuity equation 

(3.10) <9 t p + V-J = 0, 

while gauge invariance requires that the physics is unchanged by the replacement 

(3.11) A»(x) -> A»(x) - d»f{x), 

with / an arbitrary function and — (dt, — V). Equations (|3.1UI3.11|) imply 

(3.12a) q^ v = 0, 

(3.126) K" v q v = 0. 

More explicitly, one has the following relations connecting the various correlation func- 
tions: 



(3.13a) 


ljK 00 


= q j K j0 , 


(3.136) 


luK°* 


= q J K j \ 


(3.13c) 


luK " 


= q J K°\ 


(3.13a") 


uiK i0 


= q 3 K lj , 



( 3 ) The plus sign in front of the diamagnetic term is due to the fact that by using a lower index 
for Afj, the space part has a minus sign. 
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from which K° l = K i0 , K ij = K ji , and 

(3.14) u) 2 K 00 = q i K i iq j . 

The conductivity tensor may be decomposed into longitudinal and transverse components 



q l q> ( q l q 3 



(3.15) a iS = -^(tl + Uy -—j<r T 
so that ea. H3.9(l reads 

(3.16) <Ti=ihm lim ^rK 00 (q, w). 

w^o| q |^oq 2 

The charge response to a static and homogeneous external potential, e.g. the compress- 
ibility, is given as 

(3.17) ^ = _J_ l im X oo (q,0). 

o/j, e z |q|-*0 

To appreciate the physical meaning of eas. (|3. 1613. 17(1 . let us consider the phenomena-log- 
ical expression of the current for a good metal 

(3.18) J = cr L E- DX7p 

where D is the diffusion coefficient, which, under general statistical considerations, is 
related to by Einstein's relation 

(3.19) °L = e*pD. 

on 

Equation (|3.18|l may be used together with the continuity equation l|3.10|l to find an 
expression for the density-density, K 00 , response function. By taking the divergence of 
eg. (|3. 18(1 and replacing it into ea. (|3.10l) . one gets 

(3.20) (d t - DX7 2 )p = a L X7 2 qi, 

from which, after Fourier transforming, the density-density response function reads 

The above equation, of course, agrees with ea. (|3.16|) and gives the compressibility (|3.17|) 
as required by Einstein's relation. Notice that the latter, within the linear response, is 
derived from the ea. l|3.14[) . connecting the density-density and current-current response 
functions. The task of a microscopic theory, as it will be shown in the following sections, is 
to derive the expression for the current instead of phenomenologically assuming ea. <|3.18[l . 
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3'3. Response functions and Ward identities. - We begin by introducing a vertex 
function 

(3.22) A?(x, x', x") =< TtJ^ix^ix'^ix") >, 

where T t is the time-ordering operator and the average is over a statistical ensemble. In 
this subsection, for the sake of simplicity, we confine ourselves to the zero-temperature 
limit with the average taken over the ground state. We also neglect spin indices for a 
little while to keep the notation as simple as possible. When the derivative d/d x acts 
on the right-hand side of ea. (|3.22|) . one obtains two contributions. One is due to the 
derivative acting on J M and gives zero due to the continuity equation (|3.10|l . A second 
contribution comes from the time-derivative of the T t -product. As a result one gets the 
following Ward Identity: 

(3.23) -^—h^ix, x', x") = ie Six - x")G(x', x) - ie Six - x')G(x, x"), 
dx^ 

where we have introduced the single-particle Green function 

(3.24) G(x, x') = -i < T t ip{x)ip^(x') > . 

One may also consider Fourier transforms with respect to the relative coordinates x — x" 
and x'—x (The arguments of the two Green's functions in the right-hand side of ea. H3.23(l ). 
We define the Fourier transform of ea. (|3.22|l as 

(3.25) A^ss, x', x") = J dq dp <Mp-*IW-*) Jb+9/W*-*") A"(p, q) 

in terms of which the Ward Identity becomes 

(3.26) q^(p, q) =eG(p- q/2) -eG(p + q/2). 
In the above p — (e, p) and in ea. l|3.25[) 

p 

and similarly for q. The connection between the vertex function and response functions 
is obtained by introducing the truncated vertex r M defined as 

(3.27) A"(p, q) = G( P + g/2)r"(p, q)G(p - q/2). 
In terms of the response functions read 

(3.28) K^(q) = -i [dp riP, q)G( P + <?/2)r>, q)G(p - q/2), 
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where 7 M = (e, e p/m) = (e, 7*) is the bare vertex. As a check, the bare vertex is found 
by writing the Ward Identity in terms of 

(3.29) g„r"(p, ^eG^I g/2) - e CT 1 ^ - g/2), 
and using the bare Green's function expression 

(3.30) G( P ) = ' 



e-^p +i sign (|p| -p F )' 



Pf being the Fermi momentum and £ p = p 2 /2m — /i. 

We conclude this section by giving a few more consequences of the Ward identity 
(|3.26(l . First, we notice that, while in the static limit, the density-density response 
function gives the compressibility (compare ea. l|3.17[) ). in the dynamic limit we have 

(3.31) lim K 00 (0,lj) = 0. 

The above result, which is a mathematical formulation of the particle number conser- 
vation, follows from eq. l|3.26|) after taking the q-zero limit and upon integration over 
momentum p and the entire energy e range. 

Finally, by restricting the frequencies to the region (e + o>/2)(e — uj/2) < 0, and taking- 
advantage of the Ward Identity in the zero-momentum limit, one gets for the density 
response function 

KIUOM = --f /2 ^E A°(e,p;a;,0) 

J-w/2 2 ^ p 

= -ie 2 f ^ £ - (G R (P, e - u /2) - G A (p, e + lj/2)) 

J-w/2 27r p ^ 

(3.32) ie 2 ^ ( GR (P> °) ~ GA (P' °)) > 



2tt 
p 



where we made use of the fact that the sign of the frequency determines whether the 
Green's function is analytical in the upper (retarded R) or in the lower (advanced A) 
half of the complex plane as a function of frequency. By recalling the expression for the 
single-particle density of states 



N{e) = --Y j TmG R {p,e), 

TT ' 



7T 

P 

one obtains an expression for the single-particle density of states at the Fermi energy 

(3.33) MO) = N= lim -^K^(0, to), 

where the factor of 2 in the denominator is due to the spin degeneracy. 
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u(r) u ( r ') 




Fig. 6. - Self-energy in the Born approximation before and after averaging over the impurity 
distribution. The dashed line represents the average of two impurity insertions. When the 
internal Green's function line (solid line) is replaced with the dressed Green's function one 
obtains the self-consistent Born approximation. 



4. Non-interacting Disordered Electrons 

4'1. Self- consistent Born approximation. - Quite generally, non-interacting electrons 
in the presence of disorder are described by the following Hamiltonian: 



(4.1) H = J drV4(r) 



2 



V 

o^7 + u r 
2m 



where u(r) is taken as a Gaussian random variable defined by 



1 



(4.2) u(r)=0, u(r)u(r') = u 2 5(r - r') = J(r-r'). 

2t:Nqt 

In the above, N is the free single-particle density of states per spin and r is a parameter 
inversely proportional to the impurity concentration and whose physical meaning will be 
evident in a few moments. In the Born approximation 03] one has for the self-energy the 
expression (see fig El 

(4.3) ^(r,t;r',t') = ^ ~ r ° G°(r, t; r, t'). 

The superscript for G indicates that we are considering the unperturbed expression 
In Fourier space, ea. 14.3fl reads 

1 



^ sl (P' £ ) = ^r-E 



2-7riVoT *-f e-^p/ +i sign (|p'| -p F )' 

For large values of p', the real part of the sum over p' diverges, but its value does not 
depend on the energy e. This divergency is a consequence of the simple model taken 
for the scattering potential. A more realistic momentum-dependent scattering potential 
will generally cure the divergence and give rise to a finite contribution that may be 
absorbed into a redefinition of the chemical potential. The main contribution to the 
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energy dependence of the sum comes from the values of p' close to the Fermi surface. 
By following the standard procedure we pass from momentum to energy integration 

/oo />oo 
<1V- - -V„ / 'lip - 
p ~M J —oo 

where we have sent to minus infinity the lower limit of integration, since /i rs Ep is the 
biggest energy scale in the problem. Then, by residue integration, we obtain 

(4.6) £ 1 (p,e) = --Ui 5 n(e). 

To proceed in the perturbative expansion one replaces the above result into the Green's 
function and computes X 2 . At second order one has exactly the same expression as before 
except that the pole of the Green's function is now moved away from the real axis by 
the quantity l/2r. However, the residue integration does not depend on the distance of 
the pole from the real axis and one realizes that the right-hand side of eq. (|4.6|) is indeed 
a self-consistent solution which yields 

(4.7) G(p,e) 



e - C P + 27 sign (e) 



from which emerges the meaning of t as the elastic quasiparticle relaxation time. We 
stress that the consistency of the above approximation for evaluating integrals over mo- 
mentum is based on the fact that the distance from the real axis of the pole in the 
Green's function remains small compared to the Fermi energy, which corresponds to the 
condition 

(4.8) E f t > 1 

and one finds the effective expansion parameter discussed previously. The self-consistent 
Born approximation effectively selects a subset of diagrams characterized by the absence 
of crossing of impurity average lines. This sequence of independent scattering events 
leads to a ladder resummation of diagrams for the vertex part, as it will be shown in the 
next subsection. 

4'2. Vertex part and diffuson ladder. - This subsection has a twofold aim. On the 
one hand we show that the microscopic approach, at leading order in the parameter 
1/EpT, reproduces the semiclassical theory of Drude-Boltzmann. On the other hand, we 
introduce a number of technical ingredients that will be used extensively in these lectures. 
In particular, we will perform the evaluation of both the density-density and current- 
current correlation function. We begin with the density-density response function. The 
expression to evaluate reads 

/oo J, 
g- G(p + ,e + )r°(p, e ;q,c)G(p_,e_) 
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(b) 

Fig. 7. - (a) Diagram for the correlation function. The black triangle indicates the vertex part, 
(b) The vertex part is obtained as an infinite resummation of non crossing impurity lines. 

where p± = p ± q/2 and e± = e ± u>/2 and we have introduced a factor of two due to 
the spin. The truncated vertex r° is indicated by a black triangle in fig|71 The Green's 
functions appearing in the diagrams are those evaluated within the self-consistent Born 
approximation, as explained in the previous subsection. The evaluation of the vertex 
requires the evaluation of the series of ladder diagrams shown in figlHl The series of 
ladder diagrams, to be called ladder from now on, may be evaluated by solving the 
integral equation 

(4.10) Lp, P >,e(q, u>) = L<°) + J2 G (p"+, e+) G (p"_, e_) L p » p ,, e (q, «), 

p" 

where = 2jT ^ nT ■ The ladder generally depends on three momenta and two energies. 
However, as it will be shown in a few moments, the actual dependence is only on q and 
u>. This is the reason of the notation adopted. Later on we will drop the subscripts p,p' 



P +q p'+q 



+ 



+ ... = 



= L 



p p' 

8 

Fig. 8. - Ladder resummation. 
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and e. The kernel of the integral equation l|4.10|l gives a non- vanishing contribution when 
the poles of the two Green's functions lie on opposite sides of the real axis. To see this, 
let us consider the term with two impurity lines in the ladder series. It reads 

(4.11) n = ]TG(p + , e+ )G(p_,e_). 



To evaluate the integral we go to the energy variable introduced in ea. l|4.5|l . We have 
£p±q/2 = £p i P " q/2m + q 2 /8m. Since the integral is dominated by the contribution 
coming from the poles we set to pf the absolute value of p in the scalar product with q. 
The integral over the energy £ p may then be carried out by residue method as explained 
in the previous subsection. Also we notice that the condition of having poles on opposite 
sides of the real axis implies a restriction on the frequencies, i.e., e + lu/2 > and 
e — ui/2 < 0. As a result we get 

(4.12) n = d{u 2 /A-e 2 )27rN T f dQp ' 



flp 1 — ilut + ilq cos(0) 

where f2 p is the solid angle and 9 is the angle between p and q. I — vpT represents 
the mean free path due to elastic scattering. Although the angular integral may be 
evaluated exactly, in the following we will be interested in the limits lot <C 1 and (5C 1. 
which define the diffusive ( ) transport regime. It is then convenient to expand for 
small frequency and momentum the right-hand side of eq. (|4.12() and perform the angular 
integral. We get finally 

(4.13) n = 6{io 2 /A - e 2 )2nN T (l + ilot - Dq 2 r) , 

where the diffusion coefficient is given by D = vpl/d, d being the dimensionality. As 
anticipated, the ladder does not depend on the momenta p and p'. It depends on the 
difference of the incoming and outgoing momenta only and the integral equation becomes 
an algebraic one. The final solution reads 

£(0) 



(4.14) 



i - mn 

1 8(u; 2 /4-e 2 ) 
2ttN t 2 -iuo + Dq 2 ' 



This is the most important equation of this subsection. It shows how the diffusive pole 
one expects emerges from the repeated elastic scattering. In terms of the ladder the 



( 4 ) This expansion is sufficient in the low temperature regime when Tr < 1. At higher tem- 
perature with Tr > 1 one must retain the full frequency and momentum dependence of II. This 
defines the quasi-ballistic regime. A detailed discussion can be found in ref.|45|. 
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vertex reads 



e ( i + E G (p'+' e +) G (p'-' e -) L p',P, e (q^) 
p' 

fl(u; 2 /4 - e 2 ) 
— kjt + Drq 2 



(4.15) =e(l + 27riVoTL(g,w))=e 1 + 



One notices that in the zero-momentum limit the expression for T° may be obtained 
directly from the Ward Identity i|3.29fl and agrees with the above equation. With the 
vertex we may now complete the evaluation of the density-density response function. It 
is natural to split it in two parts corresponding to the two contributions in the vertex. 
The contribution, which does not contain the diffusive pole and does not have restrictions 
on the frequency range, reads 

K oo K o_o_ = _ 2ie 2 r de £ G(p+;£+) G(p _, e _) 



P 



p 

= - 2ie2 h E [G fl (P,0)-G^(p,0)] 
p 

(4.16) = -2e 2 7V , 

where the superscripts -R(A) indicate the retarded and advanced Green's functions 

(4.17) G R ^(p,e) ' 



In obtaining the above result, we have taken the static limit (frequency goes to zero first, 
and then momentum) . Given the restriction in the frequencies for the second part in the 
vertex, one could naively think that in the small frequency and momentum limit, this 
contribution would be vanishingly small. This however depends on the order the two 
limits are performed. Due to the presence of the diffusive pole in the ladder, one obtains 
for the second dynamic contribution 

K° + °_ = -2ie 2 f ^ p- Y, G( P+ ,e + ) G(p_,e_) , ] g 
J-u/2 27r p -iujT + Drq z 

(4.18) » -2e 2 N 



-iw + Dq 2 ' 



where due to the diffusive pole one may set q = and u> = in the Green's functions and 
perform the integral with residue methods. By combining together eas. (|4.16ll . (|4.18l) one 
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p 




p'+q 



P' 



Fig. 9. - Diagram contributing to the leading correction to the current-current response function 
in the dynamical region. In the ladder flows the momentum q. 

obtains the total density-density response function at leading order in Epr, in complete 
agreement with the result of ea. l|3.21[l based on the phenomenological expression of the 
current. We also note that, as expected, both the compressibility and the single-particle 
density of states are given by the Fermi gas expression Nq. 

We now turn to the evaluation of the current-current response function. For the pur- 
pose of computing the electrical conductivity, this is not strictly necessary since the Ward 
identity i|3.1bfl allows us to get a from the density-density correlation function (|3.21|) . We 
believe however that it is instructive to show how the calculation goes. The expression 
reads 



The first observation concerns the vectorial nature of the vertex. By going through the 
same steps as for the density-density response function, the integral contains a vertex 7* 
which makes the integral vanish upon angular integration. As a result the current vertex 
r* remains unrenormalized. This can be illustrated with the help of figEl Since the 
region of small q gives the dominant contribution to the integral, one can set q = in 
the Green's functions. As a consequence the two p and p' integrations in figEl decouple 
from one another and vanish for the presence of the vectorial vertex. The diagram of 
figEl which is the diffusive polar contribution to the density-density response function, 
does not contribute to the current-current response functions. However, the evaluation 
of the remaining part of the response function is more delicate as compared to the case 
of the density-density response. The reason is due to the fact that in performing the 
small- frequency limit we need to divide by oj according to ea. (|3.9|l . One then cannot 
simply take the zero-frequency limit by setting u) — before performing the integral, but 
it is necessary to make an expansion in powers of lo. After a few manipulations we get 




= -2i 




p 




|±(G*(p,0)-G^( P ,0)) 2 



(4.20) 



e 2 riQ . e 2 noT 




24 



Carlo Di Castro Roberto Raimondi 



P Pi P P2 Pj Pi+ P- P 

(a) (b) 

Fig. 10. - (a) A diagram without crossing, (b) A diagram with crossing. 



The first term, by recalling ea. (|3.5|) . cancels exactly with the diamagnetic contribution. 
The remaining term, by using ea. 1)3. 9(1 . gives the Drude formula for the electrical con- 
ductivity. 

To summarize, in this subsection we have evaluated the linear response of a disor- 
dered Fermi gas to an electromagnetic external field to the leading order in the parameter 
I /(Eft). This parameter is also the natural dimensionless coupling of the present micro- 
scopic problem which deals with the Fermi gas (characterized by the unique energy scale 
Ep) in the presence of disorder which introduces the (other) energy scale Nqu 2 ~ t _1 . 
At this order, one recovers the results of the semiclassical approach of Drude-Boltzmann. 
However, we have developed a formalism within which next-to-leading corrections may 
be investigated systematically. This will be the subject of the next subsections as far 
as non-interacting electrons are concerned. Interaction effects will be considered in the 
next section. 

4'3. Weak localization. - We have stated that the leading approximation in an ex- 
pansion in the parameter 1/ {Eft) is obtained by considering diagrams without crossing. 
We begin our discussion of the next-to-leading corrections, by showing how crossing of 
impurity lines increases the order of a diagram. A simple example is shown in fig. (|10|l . 
where both diagrams are of the same order in the impurity lines. Let us estimate these 
diagrams. By recalling the self-energy expression (|4.6|l and the Green's function l|4.7|l . 
the diagram (a) reads 



(a) = (^)W) 

1 1 

t r 



while the diagram (b) yields 



(b) 



( 1 



\2tT NqT 
1 



G(pi, e)G(p 2 , e)G(pi + p 2 p, e) 



P1.P2 



E 



(27rA r or) 2 (£ _ ^ + We)(£ _ ^ + We)(£ _£ i _ 6 _ e + 2/1 _ E £ !a + P<21±2A + iSe ) 
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Fig. 11. - Diagram contributing to the next-to- leading correction to the current-current response 
function. The Greek indices label the spin of the Green's function line. Notice that at the density 
vertex the spin is conserved. 



1 

where for brevity S e = sign(e) /2t and £i = £ pi and similarly for p 2 and p. The above 
result is obtained by noting that, due to the poles of the first two Green functions, we can 
set £i = £ 2 = in the third one. Clearly (b) is smaller of (a) by a factor l//ir. In general 
to take diagrams with crossing of impurity lines becomes a very complicated problem. 
There is, however, a subset of diagrams, the so-called maximally crossed diagrams, which 
can be evaluated. The contribution of the series of these diagrams to the current-current 
response function is shown in fjg. ll LI The corresponding expression reads 

/OO J 
^G(p,e + ) G(p,e_) L c (p,p»G(p',e+) G(p',e_), 

p,p 00 

(4.21) 

where with £ c (p, p', uj) we have indicated the series of the maximally crossed diagrams. 
As for the leading order calculation, the first step requires the evaluation of the series 
giving L c (p,p',cj). This can be done by observing that the series of the maximally 
crossed diagrams, from now on called crossed ladder, L c , may be expressed in terms of 
the direct ladder by reversing one of the electron Green's function lines, as shown in 
figlT2l This corresponds to 

(4.22) ic ( Pj p' ;CJ ) =jL (p + p ' ;CJ )__^. 



p a P p' p 



x 



pa P p' -P' (3 « 

Fig. 12. - Crossed ladder expressed in term of the direct ladder. Notice that, upon the reversal 
of the bottom Green's function line, the in and out combination of spin indices are a/3 and /3a, 
respectively. 
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Hence, the crossed ladder has a diffusive form with respect to the combination p + p' 
since one of the two Green's function lines has been reversed. This means that the 
dominant contribution comes from the region p ~ — p' and the integration over the 
two vector vertices is no longer decoupled as in the leading order case with the direct 
ladder. The direct ladder studied in the previous subsection describes the propagation 
of a particle-hole pair. The crossed ladder, when one of the lines is reversed, describes 
the propagation of a particle-particle pair and the diffusive pole is with respect to the 
total incoming momentum of the pair. This is reminiscent of the interaction scattering 
channel relevant for superconductivity and for this reason the crossed ladder is called the 
cooperon ladder. The direct ladder, on the other hand, is, in the technical jargon, called 
the diffuson. A key point to keep in mind about the difference between the diffuson and 
the cooperon is the total electric charge of the pair. While for the diffuson this is zero, it 
is twice the electron charge for the cooperon. As a consequence, in the diffusive regime, 
the cooperon is affected by the presence of a magnetic field, while the diffuson is not 
touched. We will see how a magnetic field affects the cooperon in the next subsection. 
We are now ready to complete our derivation of the correction to the current response. 
By performing the integrals (see Appendix A) in the standard way, eq. (|4.21|l becomes 



SK t] = 




2ttN t 2 -iu + £>Q 2 



= di u> — — — — 7T 

ir ^ -iu/D + Q 2 

from which and from ea. (|3.9|) . at d — 2, we get the correction to the conductivity 
(4.23) 6a = - -^r In (~) = -cr Q 2t In 

where we have resumed the physical units. The above equation represents the weak- 
localization correction and t = (A^NqDTi)^ 1 = (2irEpT is the effective small 
expansion parameter in the metallic regime. To make connection with the phcnomcnolog- 
ical scaling theory, we note that the parameter t coincides with the inverse conductance 
in units of Go divided by 2tt, t = l/(2irg). The logarithmic divergence has been cutoff at 
large Q by the inverse mean free path, i.e. the distance beyond which diffusive motion 
starts to set in. The small Q cutoff is instead provided by the inverse of some length 
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scale, L. In the zero frequency and zero temperature limit, this scale is given by the sys- 
tem size. At finite temperature, inelastic processes provide a so called dephasing length 

As discussed in the first section, the physical origin of the weak-localization correction 
is due to quantum interference. We are now in a position to appreciate the exact meaning 
of this statement. In a Feynman diagram a Green's function line describes the amplitude 
for going from one point to another. In the response functions, the two Green's function 
lines represent the operation of taking the product of one amplitude with its complex 
conjugate. Each amplitude is a sum over all possible paths so that in the product 
there will appear interference terms. The leading approximation, by restricting to ladder 
diagrams, makes an effective selection of paths. When we average over the impurity 
configurations, by connecting two impurity insertions by a dashed line, this corresponds 
to the fact that both the upper and lower Green's function lines are going through 
the same scattering center, i.e., one is considering the product of the amplitude of a 
given path with its complex conjugate. Hence, in the leading approximation there is 
no interference. When considering, on the other hand, the maximally crossed diagrams, 
one observes that the upper Green's function line goes through a sequence of scattering 
events which is exactly the opposite of the one followed by the lower line. This represents 
the interference between trajectories that are one the time reversed of the other. One 
also notices that these trajectories are made by closed loops and always come in pairs, 
due to the fact the loop may be gone around clock- or anticlockwise. 

So far we have established that the electrical conductivity acquires a logarithmic 
correction at order 1/EpT. The correction has a negative sign and signals a slowing- 
down of the electron diffusion. It is then legitimate to ask what this implies for the full 
momentum and frequency dependence of the density response function. To this end, 
we now consider the corresponding corrections with the bare density vertex 70 — e. In 
the presence of the scalar vertex, as we have seen in subsec. 4.2 already for the Drude 
approximation, the direct ladder contributes to the density response function to obtain 
the diffusive form (|3.21|) . At the order t we are now considering, many more diagrams 
have to be taken into account. Most of them cancel each other 03] an d we are left with 
those shown in fig ll3l The expression for the first diagram reads 

SK oo_ = _^ £ ST% q ^)G R (p + ) G A (p„) L c (p 1 p\u)G R (p' + )G A (p'_) ST°( q , w ). 
p,p' 

(4.24) 

In ea. H4.24[) we have already performed the sum over the frequency e, which gives rise 
to the factor ui in front. In the Green's functions, consistently with the approximations 
used up to now, we have set the frequency to zero and only the momentum argument is 
explicitly shown. The two vertex corrections <5r° at the extreme left and extreme right 
describes the direct ladders appearing in the diagram and are given by ea. (|4.15() . To 
evaluate the integrals over the momenta p and p', we note that the cooperon ladder 
gives a big contribution to the integral when p + p' is small. It is then convenient to 
introduce the variable Q = p + p' and set Q = everywhere except that in the cooperon 
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Fig. 13. - Diagram contributing to the next-to-leading correction to the density-density response 
function. 



ladder. We are then left with an integral over the cooperon ladder and a second over a 
product of four Green's functions. However, this is not yet the full story. It turns out that 
there exist further diagrams (the second and the third) that may be obtained by simply 
decorating the first diagram of fig with an extra impurity line. Such a decoration only 
adds two Green's functions and an extra summation over a fast momentum. To proceed 
we have then to integrate the Green's functions in the diagrams of fig llHI according to 
the expression 



Ih = ^G i? ( P+ )G A (- P+ )G fl (-p_)G A (p_) 



(4.25) 



G R ( P+ )G R (-p-)G A (p-)Y, G*(p' + )G A (-p> + )G R (-p>_) 
p p' 

£ G R ( P+ )G A (-p + )G A (p^ G A (-p' + )G A (p'_)G R (-p>_). 



Notice the factor 2 - for the extra impurity line. Since we are interested in the small- 
momentum limit, we make an expansion in powers of q. After a straightforward, but 
lengthy calculation, one gets 



(4.26) 



I H = AitT A N a Dc?. 



Finally, from ea. (|4. 15|) and the expression l|4.22|l for the crossed ladder, we get as a 



Disordered Electron Systems 



29 



correction to first order in t to cg. (|4.18() for K^°_ 

v y + (Dq 2 - kj) 2 7T 4^ DQ 2 - iw (Dq 2 - \uj) 2 



where 



5D 



5a D 



■2J- 



2e 2 7V ttA^o ^ £>Q 2 - iw 
Q 



agrees with the expression for 5a found in ea. l|4.23[l derived from the current-current 
response function. The above equation shows that the correction to K 00 , made of many 
different contributions, can at the end be absorbed into a renormalization of the diffusion 
coefficient Dp = D + 5D. Hence we conclude that at this order ea. (|4.18|l changes into 

(4.28) K° + °_ = -2e 2 N — U ' 



+ D R q 2 

and by virtue of ea. l|3.17|l and ea. H3.33|) the compressibility and the density of states 
are not renormalized. Only one renormalization (D — > Dr) is required in this case. 
The one-parameter scaling theory follows. Given the expression (|4.23(l for 5a, the group 
equation for the conductance g has an expansion in t of its /3-function fea. l|2.29|l % ) with 
the coefficient a = 1/ir. The critical index for the conductivity, at order e in d = 3, is 
[j, = l. The frequency cuts off the singularity in the diffusion ladder and acts in this 
transition as an external field in ordinary transitions. Its scaling index x u has the same 
value as the dimension of the Dro 2 , i.e., x u = e + 2 = cl. 

4'4. Effect of a magnetic field. - As we have pointed out in sec. 2. 3, the magnetic 
field will cut off the weak- localization corrections. To see this explictly, it is useful to 
switch to a space and time representation of the cooperon ladder, as shown in fig ll4l 
The cooperon describes the propagation of a pair of electrons that have coiciding space 
coordinates (within the spatial resolution given by the mean free path I), t, t' and rj, rf 
are center-of-mass and relative times of the electron pair. For instance, an incoming pair 
has a temporal evolution factor 

(4,29) e -i(e+^/2)(t'+>j72) e -i(e-^/2)(t'-V/2) = e ~2iet' & -iuJ V ' _ 

In this representation the cooperon reads 

9(n - n') p-|r-r'|7-D(r,-,/) 

{ ' c { ' J 2itN t 2 {2T:D{f 1 -i 1 ')) d / 2 ' 

and obeys the diffusion equation 

(4.31) ( 2 ± - DV^ LT' (r, r') = - v'W* - *')■ 
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We note that the formula for the weak-localization correction of ea. i|4.23[) involves the 
Fourier transform with respect to time of the cooperon at coinciding space points. In 
terms of (r, r'), the weak-localization correction ea. 14. 23(1 . in d — 2, is recovered as 



(4.32) 5a = -e 2 4N Q DT 2 [ * d V L^ n {r, r) = -\- In 



The space and time representation makes more transparent the physical origin of the weak 
localization correction. The cooperon propagator L2'~ v (r,r) represents the propagation 
of a pair of electrons going around the same closed trajectory in opposite directions, or 
one electron going through the time reversed trajectory of the other electron. The time 
nedeed to go around the loop is 2ij and one has to integrate over all possible values of 
r\ between r and r^. The lower limit, r, sets the time over which diffusive behavior 
develops. The upper limit, r^, sets the time over which phase coherence between the 
two electrons going around the loop is maintained. For both times, we switch to the 
corresponding lengths via the diffusion constant. 

We are now ready to consider the effect of an external magnetic field. It enters the 
diffusion equation via the minimal substitution as 

(4.33) V -> V - 2ieA(r), 

which is to be expected for the minimal substitution of the two electrons described by 
the cooperon. In fact, as the established name cooperon may suggest, this is com- 
pletely analogous to the minimal substitution that it is usually made when considering 
the Landau- Ginzburg equations for superconductivity. In the absence of the external 
magnetic field, the diffusion equation is solved via the knowledge of the eigenvalues of 
the laplacian operator exactly as for the Schrodinger equation at imaginary time. The 
analogy is made precise by saying that one may define particle-like parameters as a mass 
m* = 1/(21?) and a charge e* = 2e. A magnetic field will modify this mass and act as 
a cutoff for the singularity in a. In the presence of a uniform external magnetic field, 
the diffusion equation may be solved in terms of Landau levels. More precisely, in two 
dimensions for a magnetic field perpendicular to the plane, the cooperon at coinciding 
space points reads 

(4.34) Lr'(r,r) = g. £ e"^V» = fe^> g , ^^l, 
V ' c V ' ; 2ttN t 2 y ^ 2irN T 2 ys 1 - e -^(W)/2 ' 

n— 

where g s = e*B/ (2tt), E n = uj c (n+l/2), and ui c — 2e*DB are the degeneracy and energy 
of the effective Landau level, and the ciclotron frequency, respectively. By inserting the 
ca. (l4.34|l into the expression for the weak-localization correction of ea. (|4.32(l one gets 



(4.35) 6a(B) = -e 2 D— I drj . ' ^r = -7r^T ln 

- ' ' sinh(w c 77/2) 2n 2 h 



tanh (lOcT^/A) 
tanh (w c r/4) 
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t+ rj 12 



t'+ r|V2 



t-7i /2 



t'- n72 



Fig. 14. - The cooperon ladder in time and space representation. Notice that both particles 
in the incoming (outgoing) pair have the same (center-of-mass) coordinate r' (r). The lack of 
explicit dependence on the relative coordinate corresponds to the fact the ladder resummation 
does not depend on the fast momentum. 



In the metallic regime, where r <C r^, at low magnetic field (uj c t^, <C 1), the above 
correction reads 



(4.36) 



Sa(B) = - 



2n 2 h 



In 



4 (<&(B) 



where $o — hc/2e and $(£?) = ftL^B is the magnetic flux through a circle with radius 
the dephasing length. The effect of the magnetic field is felt for fields such that &(B) is 
of the order of the flux quantum $ - Then it suppresses the weak localization correction 
and gives rise to a negative magnetoresistance. Experimentally, by measuring the mag- 
netoconductance one may obtain the value of the dephasing time at a given temperature, 

In the opposite limit, such that uj c t^ ^> 1 (but lo c t <C 1), ea. (|4.35|l reads 



(4.37) Sa(B) = ^-ln(oj c r), 

which is no longer singular. The crossover between the two regimes occurs when the 
infrared cut off t^ 1 is replaced by uo c . This condition amounts to &(B) ~ $o> a s was 
already stressed in subsec.2.3. 

As far as the metal-insulator transition is concerned, the suppression of the weak 
localization correction corresponds to the vanishing of the coefficient a of the perturbative 
expansion of the /3-function (cf. ea. l|2.29|l ). Then, corrections arise when a crossing of two 
direct ladders is considered in current-current response function and appear at the second 
order in the expansion in power of t ~ 1/g. This, as mentioned at the end of subsection 
12 31 implies a change of the critical conductivity exponent from fj, — 1 to \x = 1/2. This 
is contrast with several experiments, as discussed at the end of this section. 
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4'5. Spin effects. - The weak-localization correction is also affected by the presence 
of magnetic impurities. These latter make the cooperon propagator massive. In the 
presence of spin-dependent scattering, it is useful to decompose the ladder into singlet 
and triplet components with respect to the total spin for the incoming pair. The latter 
is a particle-hole and particle-particle pair for the diffuson and cooperon, respectively. 
In Appendix [5] we gi ye the details on how to derive the expression for the ladder in 
the presence of spin- flip scattering by magnetic impurities. The final result, in the limit 
when the spin-flip scattering time is much larger than the elastic scattering time, t s 3> t, 
reads 

(4.38) L s 

(4.39) L T 



2itN t 2 


Dq 2 


— iuj ' 




1 




1 




2ttN t 2 


Dq 2 


— iuj -\ 


4 
3t„ 


1 




1 




2ttN t 2 


Dq 2 


— iuj -\ 


_ _2_ ' 


1 




1 




2ttN t 2 


Dq 2 


— io; H 


2 
3t 3 



(4.40) L s c = 

(4.41) L T C = 



One sees that only the singlet channel of the diffuson remains massless. While it is obvious 
that both triplet channels become massive, magnetic impurities as well as magnetic field 
break the time reversal symmetry and introduce a mass in the cooperon singlet as well. 
The a-coefficient of the /3-function expansion is again vanishing and the singlet channel 
diffuson gives rise to corrections to second order in t ~ 1/g, yielding the value [i = 1/2. 
Comparison with available experiments is postponed at the end of the section. 

Spin-orbit scattering has the global effect of reversing the sign of the quantum inter- 
ference contribution to the conductivity. In Appendix |0 we derive for both the diffuson 
and cooperon ladders the expression 



(4.42) L s 

(4.43) L T 

(4.44) L s c 

(4.45) L T C 



1 1 



2ttN t 2 Dq 2 - ico ' 

1 1 
2wN t 2 Dq 2 -iuo + ^ 

1 1 
2ttN t 2 Dq 2 - iu> ' 

1 1 
2irN T 2 Dq 2 -icu + j^-' 



We see that in this case the cooperon singlet remains massless in contrast to the magnetic 
impurities case. In the non magnetic impurity case the weak localization correction to 
the conductivity comes from both the cooperon singlet and triplet channels. As is shown 
in Appendix C, the singlet contribution is antilocalizing and amounts to 1/3 of the 
localizing triplet contribution. The latter being now suppressed, we are left with the 
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singlet antilocalizing interference contribution. In terms of the /3-function expansion, 
this means that the coefficient a changes sign. Also in this case, the comparison with the 
experiments is postponed to the next subsection. 

In all the cases discussed above, no corrections arise to the thermodynamic quantities 
like the single-particle density of states, the specific heat, and the spin susceptibility. 

4'6. A review of the experimental situation. - It is now time to compare the theoretical 
predictions obtained from the microscopic approach with the experiments. As we will see, 
there are a number of facts that suggest that a proper description of the metal-insulator 
transition cannot be obtained without taking into account the effects of the electron- 
electron interaction. 1) In semiconductor-metal alloys, as we have already pointed out, 
the metal-insulator transition is observed with a conductivity critical exponent fi = 1. 
This is in agreement with the non-interacting scaling theory. On the other hand, tunnel- 
ing measurements reveal that the density of states has strong anomalies. For instance, in 
|23|. the single particle density of states of Ge\- x Au x is measured for different values of 
the Au concentration above and below the critical concentration for the metal-insulator 
transition. On the metallic side, the density of states shows a dip at the Fermi energy, 
which in the insulating regime, at lower Au concentrations, develops in a gap. In the case 
of the Nb x Si\- x alloy, the value of the density of states, as obtained from tunneling mea- 
surements, is plotted as a function of the Nb concentration, and it is seen scaling to zero 
linearly by approaching the critical concentration for the metal-insulator transition 20 . 

2) The presence of magnetic impurities, in the non-interacting theory, should lead to a 
value (i = 1/2. However, in metal films of Cu : Mn, one finds experimentally /x = 1 |47j . 

3) In the amorphous alloy SiAu, it has been observed^HI A* = 1 both in the absence and 
presence of a magnetic field of 5 Tesla. This systems has a strong spin-orbit coupling 
which, within the single-particle scheme, should switch from an antilocalizing term in 
the absence of the magnetic field, to a transition with fi = 1/2 in the presence of the 
field. 4) Similarly, a value fj, = 1, both with and without a magnetic field, has also 
been observed in AZo.3Gao.7As : Si with a fine tuning of the electron concentration close 
to the metal-insulator transition, by using the photoconductivity effect [1^1 15(J|. Also in 
this system, it is estimated that spin-orbit scattering is relevant. 5) The experiments in 
uncompensated doped semiconductors are even more puzzling. First, as we mentioned 
in the introduction, there is the problem of the experimental determination of the value 
of the critical conductivity exponent in Si : P, whose value fi = l/2 |13llSl"lll2| has been 
questioned [HI]. In another n — type system, e.g., Si : As [22, also a value [i = 1/2 has 
been observed, whereas a close value of [i = 0.65 has been reported for a p-type system 
as Si : -B[S3], where the spin-orbit scattering is expected to be strong. The situation is 
further complicated by the experimental observation that the introduction of a magnetic 
field changes the value of the conductivity exponent to \i = 1 [SJ] as for the alloys. Besides 
the interpretation of the issues raised by the transport measurements, the experiments 
also show that there is a strong enhancement at low temperature of the electronic specific 
heatj^HI and the spin susceptibility 29 . 6) Finally, we want to comment on the problem 
of the metal-insulator transition in the two-dimensional electron gas. Although there 
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is quite a rich experimental literature on this phenomenon, very little is really under- 
stood. First, the real occurrence of a zero-temperature metal-insulator transition has 
not gathered the general consensus. In any case, from the point of view of the present 
discussion about the relevance of the electron-electron interaction in disordered systems, 
S«-MOSFET devices and semiconducting heterostructures are even a stronger case. In 
fact, if there is a metal-insulator transition, this is clearly beyond the conventional non- 
interacting scaling theory for which all states are localized in two dimensions for any 
value of the disorder. Secondly, the Coulomb interaction in these systems is expected to 
be very strong. 

By considering the effective mass of Si, m — 0.19m (mo being the bare electron 
mass), and taking a value for e ~ 11.9, at typical electron densities n ~ 10 11 cm -2 from 
ea. (|2.12|) one gets r s ~ 2 — 3, which has also led to suggest that Wigner cristallization 
may play a role |55j . A last point to make about two dimensional systems is their strong 
parallel magnetic field magnetoresistance. When the applied field gives a Zeeman energy 
g/is-B ~ Ep, the resistivity increases by more than an order of magnitude 56 . 

The interplay between disorder and electron-electron interaction will account for most 
of the questions arisen above. 

5. Interacting Disordered Electrons 

In this section, we consider the interplay of disorder and interaction|CTl EHI In 
addition to the Hamiltonian (|4.1|) we have now to consider the interaction term 



In the above summation over repeated spin indices is understood. We show that adding 
the interaction leads, in perturbation theory, to additional logarithmic corrections in 
two dimensions to both thermodynamic and transport quantities. We discuss how these 
corrections may be interpreted in terms of a disorder-renormalized Fermi liquid [SUl 1611 
1621 1631 Ifjl) . To achieve such a goal, a key step is the identification of which additional 
parameters besides t are required to take into account the interaction in the scaling de- 
scription of disordered systems In the next subsection, we will see how the above 
mentioned logarithmic corrections arise in the single-particle density of states, the electri- 
cal conductivity (which was the only quantity affected by disorder in the non-interacting 
case), the specific heat, and the spin susceptibility. The last two will be evaluated via 
the correction to the thermodynamic potential. We will recognize how disorder selects 
particular regions of transferred momentum in electron-electron interaction giving rise to 
the effective couplings of the theory. In the following subsection, these will turn out to be 
related to the Landau scattering amplitudes. Once the skeleton structure of the theory 
is developed, via the Ward identities we will identify the singular scale-dependent terms 



(5.1) 




p,p',q 
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(d) (e) (f) 

Fig. 15. - Diagrams for the correction to the Green's function to the lowest order in the inter- 
action, (a) and (c) are exchange-type diagrams, while (b) is Hartree-type. Interaction is shown 
as a thick dashed line. The effective interaction diagrams ((d),(e),(f)) are obtained by cutting 
the internal Green's function with momentum p q. The impurity "dressing" of the basic 
interaction diagrams selects three different regions of transferred momentum. In diagram (a), 
the momentum q ~ 0, that flows in the ladder, is also present in the interaction. (See (d)). In 
diagram (b), instead, the different position of the impurity ladder makes the momentum flowing 
in the interaction line, p' — p" , unrelated to q flowing in the ladder. (See (e)). Finally, the 
crossed ladders in diagram (c) select, in the small-g limit, the interaction in the Cooper channel. 
(See (f)). 



with the Landau parameters themselves and derive the corresponding group equations 
in d = 2 + e. 

5*1. Perturbation theory and the search for the effective couplings. - 

5*1.1. Density of states. We begin by considering the diagrams for the Green's 
function to lowest order in the interaction. These are shown in fig ll5l The basic diagrams 
are the usual exchange and Hartree contributions. Here, these two basic diagrams are 
"dressed" by the presence of disorder. As before, this is done by averaging over the 
impurity configurations. To illustrate how the mechanism now works, let us consider 
first the exchange diagram (a) in fig ll5l Its expression, before the impurity averaging, 
reads 

6G{r,r';e n ) = -T^2 / dr i dr 2 G(r, ru e n )Q(ri, r 2 ; e„ - oj m ) 
(5.2) x £(r 2 ,r';e n )V(ri,r 2 ;w m ) 



36 



Carlo Di Castro Roberto Raimondi 



where we use the Matsubara formalism. e„ = (2n + 1)ttT and uj m = 2to7tT are fermionic 
and bosonic frequencies. The above expression is written in real space and is valid for 
a given impurity configuration. The electron-electron interaction is assumed to depend 
on frequency, since this will allow us to include retardation effects usually introduced 
by screening effects. As a matter of convention, we use calligraphic letters to indicate 
quantities depending on the Matsubara frequency. As a general strategy, we first per- 
form the analytical continuation to real frequencies and then average over the impurity 
configurations by exploiting the technique developed in the previous section. Upon ana- 
lytic continuation, ie„ — ► e + iO + 02]j from ea. (|5.2|) the correction to the retarded Green's 
function becomes 

5G R (v,v';e) = - J ^ J drj J dr 2 G R (r,r r ,e) [b(uj)G R (r u r 2 ; e - to) 

{V R (r 1 ,r 2 ;u;)-V A (r 1 ,r 2 ;u J )) 

(5.3) + f(u - e)V R {n, r 2 ; cj)(G R (r u r 2 ; e - w) - G A (r u r 2 ; e - u))] G R (r 2 ,r'; e) 

where b and / are the Bose and Fermi function, respectively. V R ' A are the analytical 
continuations of the dynamical interaction V. The average over the impurity config- 
urations has two effects. First, each Green's function is replaced by expression (|4.7|l 
obtained within the self-consistent Born approximation. Secondly, one has to insert di- 
rect and crossed ladders wherever possible. This gives the leading approximation. We 
have learned that the insertion of ladders is only possible when Green's functions have 
poles on opposite sides of the real axis, i.e., when the ladder connects a retarded and an 
advanced Green's function. As a result we need to keep only the last term of ea. l|5.3ll . 
The correction to the density of the states then reads 

SN(e) = --lmJ2 / |4 /(c- e )^(q,c)L 2 (q, W )G i? (p,e)G' 4 (p^q,e- W )G«(p,e) 

q,p ~°° 

(5.4) x £ G R (p', e)G A (p' - q, e - u) £ G R (p", e)G A (p" - q, e - w) 

p' p" 

where the average over the impurities has restored the translational invariance and we 
have gone to the momentum representation. The q-integral is dominated by the diffusive 
pole of the impurity ladder. This implies that in the remaining momentum integrals one 
can perform, as usual, a small-q expansion. As it was also remarked in the non-interacting 
microscopic theory, we will perform the fast momenta (flowing in the Green's functions) 
integrals first, since they contribute only to the coefficient of the singular terms arising 
from the integration over the slow momenta flowing in the ladders. We then set q = 
in all the p-, p'-, and p"-integrals, which can be carried out in the standard way with 
the residue method and, using the results of Appendix one gets 

(5.5) SN(e) = ixmV f° d^ f(w - e) ^ 



Disordered Electron Systems 



37 



At zero temperature, the Fermi function becomes a step function giving the condition 
u> < e. In the case of short-range interaction V R remains finite in the small frequency 
and momentum limit. Due to the presence of a double diffusive pole, the integration in 
the region Dq 2 r, \lo\t < 1 over frequency and momenta gives a logarithmic divergence 

(5.6) = N V R (0,0)thx\er\ = ^ln|er|. 

The Hartree diagram (b) may be evaluated in a similar way. Its expression, after impurity 
averaging, reads 

SN(e) = -ImV2^ 2 / f(u - e)i 2 (q, u) 

* J-™ 27T! 

(5.7) x G B (p, e)G A (p - q, e - u)G R (p, e) 

where the relative minus sign and the factor of 2 are due to the extra fermionic loop in 
the Hartree diagram. The interaction parameter V% takes into account the scattering at 
large angle across the Fermi surface, as shown in diagrams (b) and (e) of fig 1 151 



±- V G R (p',e)G A (p'-q,e-u;)V R (p'-p",0) 

v o ,,, 
p ,p 

(5.8) x G fl (p",e)G A (p"-q,e-c). 



v-2 = 

N „ 
p ,p 



We then see that the presence of the diffusive pole of the ladder, by making the small 
q region more relevant, effectively selects the electron-electron scattering at small, V\ 
(exchange contribution), and large momentum transfer, V2 (Hartree contribution). A 
similar analysis can be carried out for the exchange contribution with crossed ladders 
(diagram (c)) and one has a third parameter V3, 

v * = jt ^ gR{p '' e)GA(q " p '' 6 " w) ^ (p ' + p "' 0) 
p',p" 

(5.9) x G fl (pV)G A (q-p", e -^). 

Actually, there is also the Hartree contribution with crossed ladders, but it contributes 
with minus twice the same scattering amplitude. The total correction to the density of 
states reads then 

(5.10) ^M = (y 1 _2V r 2-V 3 )t]n|er|. 

We notice that V3 corresponds to the interaction in the Cooper scattering channel(see 
fig- 030 ■ Since its presence does not change the results qualitatively, to simplify the 
exposition, we let it drop from our subsequent discussion. On a formal level, one may 
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assume the presence of a small magnetic field which, as we have seen, by introducing 
a mass, kills the singularity in the Cooper channel. The same selection of relevant 
momenta appear also in the perturbative calculation of the electrical conductivity and 
thermodyamic potential, as we are going to show. Before leaving the density of states, 
we notice that the above results are modified in the presence of long-range Coulomb 
interaction, which leads to log-square singularity in two dimensions. Details are provided 
in Appendix IDl 

5T.2. Electrical conductivity. The impurity-averaged diagrams responsible for the 
corrections to the conductivity are obtained in Appendix [E| The procedure is similar 
to that followed in the case of the density of states, but there are many more diagrams. 
For this reason, the detailed derivation of how to perform the impurity average and the 
integration over the fast momenta is left to the Appendix [E| Here we give directly the 
final result for the exchange diagram containing the interaction amplitude V\ before the 
last integration over the slow momentum and frequency. We have 

( 5.n) fa _-£2 El tf 



y*(q,^) 
(Dq 2 - iw) 3 



where the dimensionality factor d comes from the angular integration over q. In two 
dimensions, in the case of short-range interaction, ea. H5.11f) yields |5lfll59| 

(5.12) ^-^- 2 ^H^)' 

where, as for the density of states, the factor — 2V2 takes into account the contribution 
coming from the Hartree diagram. In the case of Coulomb long-range forces, , at small 
q, the screened interaction introduces an extra singularity (see Appendix[D| V R (q,ui) ~ 
— 2irie 2 ui/ (Dq 2 ). In contrast to the case of the single-particle density of states, the 1/q 2 
singularity arising from the interaction is compensated from the additional q 2 factor in 
the integrand. By using ea. (|L).l|l for V R (q,uj) one obtains for the exchange contribution 

^ fo = -2^ ln (r7)< 

which has the same form as the weak localization correction. The same will hold for 
both short- and long-range case in the renormalized perturbation theory as we shall see 
in the next subsection. 

5T.3. Thermodynamic potential. Contrary to the non-interacting case, the interplay 
between disorder and interaction introduces singular corrections to the specific heat and 
spin susceptibility. We present here the derivation of the correction to the thermodynamic 
potential. The effective diagrams are shown in fig. We have to carry out the 
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a 



(a) (b) 

Fig. 16. - Diagrams for the thermodynamic potential, (a) Hartree. (b) Exchange. The ladder 
arising from the impurity average is already shown. Notice that in the Hartree diagram the spin 
on the two lines of the ladder may differ from zero, whereas it is always zero for the exchange 
one. 



integration over the fast momenta flowing in the Green's functions. After doing that, we 
obtain 



(5.14) 



Att = -(V 1 -2V 2 )T^2J2 



U) m q 



Dq 2 + \ur, 



In Appendix[Fl we provide the details of the derivation. The Matsubara frequency sum is 
limited by |w m |r < 1. To relax the constraint in the sum, we introduce a cutoff function 



and the sum runs between minus and plus infinity. We may then perform analytical 
continuation to get 



1 



(5.15) Ail = --(V 1 -2V 2 )J2 



duj , 



-b(oS)Im 



-2iw 



— icj + t 1 J Dq 2 — ilu 



In two dimensions, the sum over q may be done and ea. (|5.15|) becomes 



An = - 



4n 2 D 



dto ui b{u>) ln(|a;|r) 



(5.16) 



-t(Vi - 2V 2 ) ^ — ln(Tr), 



where we have dropped non singular terms in temperature. The thermodynamic potential 
then acquires a logarithmic correction that implies for the specific heat|66| 

(5.17) SC V = C Vfi t{Vi - 2V 2 ) ln(Tr), 

where Cv.o = (2tt 2 NqT)/3 is the non-interacting value. 
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In order to evaluate the spin susceptibility, we must include the Zeeman coupling. 
As shown in the eq. (|G.2j) of Appendix only the ladder corresponding to total spin 
±1 of the incoming particle- hole pair are affected by the magnetic field via the Zeeman 
energy uj s — gfj,gB. In this respect, we notice that only the Hartree diagram of fig llBI 
contributes, since in this case the total spin of the particle-hole pair is given by the 
combination a — f3, with both indices running over ±1/2. In the exchange diagram 
the total spin of the particle-hole pair entering the ladder is a — a and is always zero. 
The combination of the two diagrams may be arranged as the sum of a singlet and 
triplet component with respect to the total spin of the particle-hole ladder. Hence, the 
magnetic field only affects the triplet component with value M = ±1. For the purpose of 
isolating the magnetic-field dependent contribution to the thermodynamic potential, it 
is convenient to write the difference, with and without the magnetic field, of the triplet 
as 



AO, 



ui m Mil q 



i 



W m q 



Dq 2 + \oj m \ — \Muj s sgn(uj) Dq 2 + \uj„ 
2{Dq 2 + \u m \) 2 



{Dq 2 + \u m \) 2 +LU 2 Dq 2 + \LU ri 
-7TtN V 2 Tj2Wm\ln 

TTtN Q V 2 LU 2 s T V — 

* J C,).__ 



w m >0 



(5.18) 



—z—oj s ln(Tr) 



Now a few words concerning the steps leading to the final expression of ea. l|5.18|l . Since 
the factor \tu m \ in the sum excludes the term with oj m = and the expression in the 
sum is even in u) m , we have rewritten the sum as twice the sum over the strictly pos- 
itive frequencies. Then we observe that the smallest frequency is 2ttT, which allows 
us to make a small magnetic field expansion u s < T. From ea. l|5.18|l one finally gets, 
by differentianting twice with respect to the magnetic field, the correction to the spin 
susceptibility 



(5.19) x=-Xo2^ 2 ln(Tr), 

where we have introduced the non-interacting value \o — 2A r o(5/is/2) 2 . No correction 
is instead found for the compressibility. 

The two parameters V\ and V 2 that appear in the perturbative expression of the 
density of states (|5.1U|) , conductivity (|5.12|) , specific heat l|5.17|l , and spin susceptibility 
(|5.19() are the natural candidates for the additional running couplings to be used with 
the dimensionless resistance t to obtain the renormalized perturbation theory in 2 + e 
dimensions. 
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5'2. The renormalized perturbation theory and effective Fermi-liquid description. - In 
this subsection we show how the perturbative results derived in the previous one may be 
generalized to all orders in the interaction and how the effective scattering amplitudes 
arc related to the Landau parameters. We then develop the renormalized perturbation 
theory for the various response functions by making use of the Ward identities that 
implement the conservation laws for charge, spin, and energy. Finally, we go back to the 
perturbative results which will be generalized to all orders in the interaction and derive 
the group equations. 

5'2.1. Effective scattering amplitudes and Landau parameters. According to the 
discussion of the previous section, the relevant interaction terms in the Hamiltonian are 
reduced to 



where the primed q-summation will be confined to small q values as implied when per- 
forming the disordere averaging. 

Up to now our discussion has been limited to the first order in the interaction. How- 
ever, we do not want to confine our theory to small interaction and the only true ex- 
pansion parameter must be the dimensionless resistance t = l/(2i:g). The good metal 
condition, g ^> 1 ( Epr/h 3> 1), implies that the disorder only affects electron states 
within a small distance Ti/t away from the Fermi surface. Under these circumstances, 
one may go beyond the first-order interaction correction by replacing V\ and V% with the 
Fermi-liquid scattering amplitudes Ti and r 2 , whose lowest order diagrams are depicted 
in fig(T7| 

We note that in the absence of spin-flip mechanisms, the total spin of two colliding 
particles is a conserved quantity. It is then convenient to introduce the singlet T s and 
triplet r t scattering amplitudes. After the selection of the relevant momentum transfer 
terms (cf. eas. l|5.Gll - <|5.8[l ') the corresponding spin structures can be decomposed as 

(5.21) — r 2 <S Q /3(5 7 5 = - -r 2 )(5 a ,5<5/3 7 — -T 2 cr a s ■ CX3 7 

where we have used the identity 



(5.20) F, = EE[ 



^l a ap a /3p'+q a /3p' a ap+q 



t f 

+ ^aLpOgpZ+qa^p+qaap' , 



p,p' q 



(5.22) 



We then define 



(5.23) 
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+ q / ~ \ p'+q 
p V7 p' 



p p 



Fig. 17. - Lowest order diagrams for the small (Ti) and large (F2 ) scattering angle. Notice the 
different momenta flowing through the interaction propagator. Also the two diagrams have a 
different spin structure. 



The scattering amplitudes r s and T t are related to the Landau Fermi-liquid parameters 
F° and F a °by[lU 

IF IF 
(5-24) T s = s —, T t = -- 



2A^o 1 + F° ' 27V 1 + F° ' 

From now on, when necessary, JVq is assumed to include the Landau effective-mass correc- 
tion. In terms of the Landau parameters, compressibility, spin susceptibility and specific 
heat are given by 

Tip = IT^° = 2No{1 - 2NaTs) 55 2NaZ " 
x = = xo(i + 2iv r t ) = xoz? 

(5.25) Cy = Cy i0 = CvfiZ , 

where Z° will be different from one in the presence of disorder and is here introduced for 
symmetry in the equations. 

5'2.2. Renormalized response functions. We begin our discussion of the interplay 
between interaction and disorder in the renormalization of the response functions, by 
considering the density-density response function. As in the case of the non-interacting 
theory of eas. (|4.1bj) . l|4.18[) . we split the response function in static and dynamic contribu- 
tions. The effect of the interaction may be understood in terms of a skeleton perturbation 
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+ 





+ 



(a) 




(b) 



+ ... 



(c) 



Fig. 18. - Skeleton structure of the perturbative evaluation of the response function in the 
presence of interaction. The plus and minus signs indicate the sign of the energy flowing in 
the line. We remind that the ladder insertion is only possible between a pair of retarded and 
advanced Green's functions. The first diagram (a) is the one used in the non-interacting theory. 
Diagram (b) represents how the interaction "dresses" (a). It leads to i) a "dressing" of the vertex, 
here indicated as a black triangle, ii) renormalization of the ladder. Interaction appears explicitly 
in diagram (c), which has to be considered together with all the other diagrams, indicated by 
the dots, obtained by the infinite resummation of the interaction. Depending on whether we 
consider the charge or spin response function one has the vertex A s ,t and interaction T S: t- For 
the energy response function, there is also a vertex Ae, but there is no infinite resummation of 
the interaction for the reasons explained in the text. 



theory, as shown in fjgll8l The first diagram (a), due to the ladder insertion, gives the 
dynamic contribution in the non-interacting case (cf. ea. H4.18fl '). Its "dressing" due to 
the interaction is represented by the diagram (b), whose contribution can be written in 
the form 



which generalizes ea. l|4.18[i . We now discuss it in detail. We begin with the ladder, which 
is the most important ingredient of the theory. The ladder in fact is the effective prop- 
agator of the diffusive mode, responsible for the singularities appearing in perturbation 
theory. It requires a wave function renormalization a frequency (effective external 
field) renormalization Z, and a renormalization of the diffusion constant( 5 ) 



(5.26) 



icj2e 2 ^oC 2 A^ 
Dq 2 - \ujZ 




SI 



(5.27) 



11 1 c 2 



2tt N t 2 Dq 2 - iuj 2ttN t 2 Dq 2 - iZoj ' 



( 5 ) We leave the same symbol as before for the interaction-renormalized diffusion constant to 
keep the notation simple. Everywhere in this subsection this is understood whenever we are 
dealing with the renormalized ladder. 
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Fig. 19. - Dynamical resummation of the interaction. 



In Appendix [H] we show that the interaction corrections to the ladder do not destroy 
its diffusive pole behavior and indeed lead to the renormalized form of ea. l|5.27[) . Fur- 
thermore, the logarithmically singular corrections to the single-particle density of states, 
specific heat, and conductivity which appear in two dimensions can be absorbed in the 
above three renormalizations £, Z, and D, respectively. Here, by exploiting the con- 
straints given by the general conservation laws embodied in the Ward identities, we will 
be able to directly express the ladder renormalization parameters in terms of physical 
quantities. 

The vertex A s is one in the non-interacting case. In the interacting case, it repre- 
sents the vertex which, when multiplied by A^°_, gives the total dynamic part of if 00 , 

which includes also terms ending with two advanced (++) or two retarded ( ) Green's 

functions. The vertex A a is irreducible for cutting a ladder propagator. 

Besides "dressing" the non- interacting diagram, interaction leads to new diagrams as 
diagram (c) of fig l!8l which gives 

00 _ 2 uj2N ( 2 A s ^ Lu2N ( 2 A s 



(5-28) K% = e' ~,\- T . 

y ' Dq z — vjjZ Dq z — iojZ 

By keeping in mind that the order in the expansion parameter t is determined by the 
number of integrations over the momenta flowing in the ladder propagator, we can re- 
place, without changing the order in t, the scattering amplitude r s by its screened form 

(5.29) T s (q, u) = T s - T s i^(2^V r) 2 L(q, cj)T s (q, cj) = T s ®f ~ ^ , 

Zir Dq z — iZ s oj 

obtained by an IIP A- like infinite resummation (shown in fig ll9f) and using ea. l|5.27|l . 
In ea. (|5.29|l Z s = Z — 2( t 2 NoT s will turn to be the expression dressed by the disorder 
of Fermi-liquid renormalization of the compressibility Z^. If we insert the dynamical 
amplitude (|5.29|) in the expression (|5.28l) for the first interaction correction and combine 
it with the contribution (|5.26|) of the first diagram, we arrive at the final form of the 
response function which resums all the infinite series of diagrams indicated by dots in 

fig- El 

(5.30) AT = _ ^oC 2 A S ; ^ 00 = _ g2 an + 

Dq z — iu}Z s an 
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We are now ready to make use of the Ward identities (|3.31|) . l|3.33[) . where now K+L is 
given by ea. (|5.30(l . and we get 

(5.31) N C ^ 



N Z s ' 
dn C 2 A 2 

(5.32) p = 2N ^. 

dfi Z s 

We can now follow two alternative routes leading to the same result. We can first 
remember from Appendix El that N/Nq coincides with £. It then follows from ea. (|5.31(l 
that (A s = Z s , which when used in ea. (|5.32|) gives 

(5.33) — = 2N Z S . 

Alternatively, it can be shown that the compressibility has no logarithmic corrections, 
i.e., is given by the Fermi- liquid value in the absence of disorder with Z s = Z® = (A s . 
This means that, although both A s and £ have logarithmic corrections, the combination 
(A s does not. The final expression for the density response function has the conserving 
form 

(5.34) K^{ q ,u) = -e 2dn °^ 



d\i D c q 2 



where the charge diffusion constant D c = D/Z s has been introduced. 

A similar analysis can be done for the spin susceptibility and the specific heat. First 
we note that the dynamical scattering amplitude in the triplet channel reads 

(5.35) T t (q,u) = T t +T t i^(27:N a T) 2 L(q,cj)T t (q,u;) = T t - ' ^ , 

Z7r Vq z — lZtui 

where Z t = Z + 2No( 2 T t is the expression dressed by the disorder of the Fermi-liquid 
renormalization of the spin susceptibility Z® = l+2Aor f , which as for the compressibility 
is given by the static limit (u) = 0, q — * 0) of the spin-spin response function x(q, to). 
To analyze the spin-spin response function, one introduces, in analogy with the density- 
density response function, the i-th component of the spin density 

(5.36) S\x) = 

and the associated spin-current density S\ By proceeding as in the charge density case 
of ea. l|3.22|l . we introduce a vertex function 



(5.37) A^(x,x',x") =< T t S l ^{x)^ a {x')^l{x") 



>, 
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where the first upper index indicates which component of the spin density we are dealing 
with while the second one, fj,, distinguishes between density and current component. In 
the presence of a magnetic field, one may derive for A l,M the following Ward identity 

(5.38>fcA$(p, g) = ^ [a% G (i {y - g/2) -G a (p + g/2) a^} + i ( , /:1 ^.V/>. g), 

where etjk is the full antisymmetric tensor and we have assumed the magnetic field in the 
^-direction. We now consider the consequences of eg. (|5. 38[) in the vanishing momentum 
limit. Instead of working with the x and y components, it is convenient to switch to the 
circularly polarized ones defined as 

(5.39) A^° = i^'O + A 2 ' ), 

(5.40) A iT '° = ^ (A 1 ' - A 2 ' ) , 

corresponding to the spin-density in the M = =pl triplet channels, respectively. By 
including the M = channel of the triplet corresponding to A 3:0 , one gets the Ward 
identity l|5.38|) in the form 

(5.41) ( W + M Wa )A$ (e, u) = ^ [a% G,(e - u/2) -G a (e + w /2) <r^] , 

where <r M are defined as in eq. (|5.39|) , (|5.4U|t and we have dropped the explicit dependence 
on momentum p both in the vertex and in the Green's functions. ( 6 ) The dynamical 
resummation of the skeleton structure, analog to ea. (|5.3U|l for if 00 , with T s replaced by 
T t and dn/dfi or 2N by x, gives 

(5.42) X l I _(q,u 1 )= '^ X f kt M7 , X M (q^)^ X + xl I _(Q,u J )A u 

Dq z — iuZt — lM Zh IjJ s 

where we have introduced a renormalization factor for the Zeeman energy, Zh- Equation 
(|5.41l) together with the analog of eq. I|3.31fl for the spin-spin response function in the M- 
th channel leads to 

(5-43) 3^(0,0,) = x ^^_, 

U! + MbJ s 

(5.44) x¥-(0,") = -Xo N ' 



N oj + Mlu s 



Indeed, in the limit of zero magnetic field ea. (|5.43|) gives the total spin conservation 
by the vanishing of the response function at finite u> as q goes to zero, while ea. (|5.44() 
reproduces the single-particle density of states in term of the dynamical part of the spin 



( 6 ) Notice that the Green's function gets a spin label in the presence of Zeeman coupling. 
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response function. Finally, by making use of the Ward identities (|5.43ll . (|5.44() in l|5.42[l . 

one obtains 

(5.45) X M (g,c) = X n D 2 Sq2 -7 n 1 /r MUJ l s ,Zt = CA*,x = X oZt,Z H = Z u D s = D/Z t . 

Notice that the result Zh = Z t 68 implements, to all orders in the expansion in the 
parameter t, the Fermi-liquid zero-order result. 

The energy-energy response function \e can also be decomposed according to the 
eq. (|o"4*2T) . where \ is replaced by \e and Z t and 2N K 2 t by Z and CVnTA| |6lfl 170] . ( 7 ) 
The analogous Ward identity gives 

(5.46) Z = (A E , C v = C v ,oZ, D E = D/Z 

and the frequency renormalization Z of the ladder is identified with the specific-heat 
rcnormalization. In the presence of Coulomb long-range forces, to avoid double counting, 
one has to subtract the statically screened long-range Coulomb Tq from the full singlet 
scattering amplitude entering the ladder resummation for the density-density response 
function. Hence, T s — > T s — Tq in Z s , where 

(5-47) r (g,w = 0) = ita/\o /£i "?^o 



1 + V c (q)dn/dn H ' dn/dfi 
As a consequence of ea. l|5.32|l . from Z s = Z — 2NqC 2 (T s — Tq) we derive the constraint 
(5.48) Z = 2N ( 2 T S . 

In this way we have completed the general formulation of the effective renormalizcd 
Fermi-liquid theory. 

5 '2. 3. Derivation of the group equations. We now come back to the perturbative 
expressions of the electrical conductivity (|5.12|) . specific heat (|5.17|) and spin suscepti- 
bility H5.19fl . These expressions were derived to first order in the interaction. However, 
we have seen that we may relax this condition in two ways. First, we may use the Fermi 
liquid scattering amplitudes. Second, we can make an infinite dynamical resummation 



( 7 ) In contrast to the density and spin response functions, in the energy response function the 
renormalization parameter does not require additional terms due to the interaction besides Z, 
i.e., the Fermi-liquid renormalization, the analog of diagram (c) of fig. 118H . is missing. This 
occurs since the interaction separates the integration at the two vertices in the response function 
diagrams. Due to the presence of the energy in the thermal vertex, each integration over the 
energy contributes at least a term with the second power of the frequency. Hence the terms due 
to the dynamical resummation of the interaction give rise to negligible contributions going with 
the fourth power in the frequency. 
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(eqs. H5.29j l. 15.35j l). For the electrical conductivity, by inserting the dynamical scattering 
amplitudes in ea. (j5.11jl one has 



5a 
a 



(5.49) 



2 
2tti 

1 + 



3 1- 



00 j 9 ( u / n 

aid— to coth — 



Z a Zo 

-7i — In — 

2iV C 2 r s Z 

Zt Zt 

2N er t n z 



E 



Dq 2 (T s (q,m)~3T t (q,L0)) 
(Dq 2 - iZuj) 3 



In (Tr). 



which coincides with SD / D found from the ladder renormalization in Appendix[HJ In the 
case of Coulomb interaction, when 2N ( 2 T S = Z, the singlet part of ea. (|5.49|) reproduces 
ea. (|5.13j) . Furthermore for small T s and T t (Z = 1) one recovers the first order short 
range interaction result of ea. (j5.12j) . However, in the short-range case, by allowing the 
group equations to flow, Z s = Z® is invariant, while, as we shall see, Z and 2Nq(^ 2 T s 
diverge and the singlet strength becomes again universal in Ea H5.49j) and equal to one. 

The full expression for the correction to the thermodynamic potential is obtained by 
introducing in cas. lj5.14ji . H5. 181) the singlet and triplet dynamical scattering amplitudes 
and the rcnormalized ladder, including the magnetic-field renormalization Zh = Z t . For 
the magnetic-field independent part one then gets 



^E 



quJn 



dA 



iv c 2 r s |w m | 



Dq 2 + (Z - A27VoC 2 r s )|w„ 



(5.50) 



37V C 2 r t |^, 



Dq 2 + (Z + \2N ( 2 T t )\u Jm \_ 
while the field-dependent contribution reads 

N C 2 TM 



(5.51) 



T E / dA E 



Mil 



Dq 2 + (Z + \2No( 2 T t )\uj m \ - iM(Z + 2iV C 2 r t )w s sgn(w) 



Dq 2 + (Z + \2N ( 2 r t )\u; n 



In the above equations, due to the presence of the dynamical resummation of the in- 
teraction we used the standard trick [44j of multiplying the interaction by a parameter 
< A < 1. However, this must not be introduced in the amplitude present in the mag- 
netic field insertion since it will generate spurious diagrams. As a result, the corrections 
to the specific heat and to the spin susceptibility are 



(5.52) 
(5.53) 



SC V = C v , t(N C 2 T s - 37V C 2 r t ) ln(TV), 



5x = -Xo4tAT C 2 r t Jln(Tr). 
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The renormalization of the specific heat may be interpreted as the renormalization of 
the quasi-particle density of states Nqp = ZN$. Accordingly, the renormalization of the 
spin susceptibility must contain both the renormalization of the quasi-particle mass and 
of the Landau parameter F®. To show this we write 

(5.54) X = XoZ^ = XoZ (l + 2N ^ Tt ^j = XoZ (1 + 7t ) 

with 7 t = 2N ( 2 T t /Z being the renormalised Landau static amplitude. We note that 
£ 2 is always associated either with T s or T t and drops out from the following group 
equations. It however renormalizes the single-particle density of states, which, in the 
interacting case, becomes scale dependent even though in a complicated way. Let us 
define the flow variable s = — In Tt so that s — > oo corresponds to the infrared limit. 
Then we have 

dZ t „, 

(5-55) — = - -Z(l-37t), 



(5.56) ^ = 2tZ 7t (l+ 7t ). 



According to ea. (|5.54|) one has 



ds 2' 
dZ* 
ds 



(5-57) — = Z— + (1 +lt)-r- , 

as as as 

from which, by using eas. l5.55|) . H5.56fl . one obtains 
(5.58) £-50 +■»>"■ 

in complete agreement with the explicit diagrammatic evaluation of the disorder induced 
corrections to the scattering amplitudes lr!2"] . By writing the correction to the con- 
ductivity in ea. (|5.49|) in terms of 74, the dependence on Z drops out. One gets 



dt o 

5.59) — = t 2 

as 



3(l-i^ln(l + 7 t) 

It 



By resuming the weak-localization contribution, one obtains a term identical to the sin- 
glet contribution (the first one in the square brakets). The two terms although identical, 
have therefore a complete different origin and this shows up in the presence of a magnetic 
field which kills the weak-localization contribution and does not affect the singlet one. 
Equations (|5.58(l . l|5.59|l together with ea. (|5.55|l are the renormalization group equations 
at one-loop order for the problem of interacting disordered systems at d = 2. Their 
analysis is the task for the next section. 
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6. The Renormalization Group equations 

By resuming our maritime metaphor, we finally have land in sight. In the present 
section, we approach the end of our trip by examining the consequences of the physical 
picture we have developed in the previous sections and briefly compare [2] the results 
with the experiments. To this end we discuss here in some detail the solution of the 
renormalization group (RG) equations for the inverse conductance t, triplet scattering 
amplitude 74 and the parameter Z. We begin our discussion with the general case when 
there is no magnetic coupling in the system. In this case the RG equations read 



l + 3(l_i±T*l n (l +7t ) 

It 



(6.2) g-W, 

(6.3) ^ = -| Z(1 -3 7 «), 

where in eq. (|6.1|l we have added the contribution due to the bare dimension due to Ohm's 
law, e = d — 2. The first observation is that eas. (|6.1|) . H6.2|) do not depend on Z. After 
solving for t and 74 one may successively solve ea. (|6.3[l for Z. 

Let us consider first the case d — 2, i.e., e = 0. Equation l|6.2|l for -f t implies a 
continuous growth. By integrating it between sq and s, one has 

1 1 1 r s 

(6-4) ; rr = - n / ds'i(s'), 



from which one sees that jt diverges at a finite value, s c , of the flow parameter: 
(6.5) l = ±(l+ 7t (s Q )) J\s%s'). 

The ea. l|6.1[l for t says that after an initial increase for not too large 74 (so), the growth 
of 74 makes the triplet contribution, which is antilocalizing, the dominating one. As 
a result, t goes through a maximum. In fig l20l we show the RG flow in terms of the 
variable t/(l + t) and 7t/(l + 7t). For all the RG trajectories 74 = 00 at some finite 
value s c , which depends on the initial values. Due to this, one cannot seriously trust 
the above equations quantitatively. Nevertheless, the physical indication of some type of 
ferromagnetic instability is rather clear due to the diverging spin susceptibility associated 
with 7 t . The appearance of a finite length scale may indicate a formation of local magnetic 
moments on the same scale. Furthermore, the dominating antilocalizing effect of the 
triplet while t remains finite strongly supports the possibility of a metallic phase at low 
temperature |731 172*1 1671 l?4l I75| . in contrast with the non- interacting theory based on WL 
only. Indeed, this metallic phase in d = 2 has recently been observed (see refs. in |34| ) . 
In any case, both experimentally and theoretically, it is not clear whether a possible 
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ferromagnetic phase occurs before the transition to the insulating phase. Almost all the 
experimental information is based on transport measurements and the spin susceptibility 
is obtained indirectly. Very recently, a new method for measuring directly the spin 
susceptibility in a two dimensional electron gas has been invented and the first result 
suggests that, although there is a spin susceptibility enhancement, no ferromagnetic 
instability is observed [75], 




0.4 0.6 
t/(1+t) 

Fig. 20. - The RG flow, for d — 2, in terms of the variable t/(l + t) (a;-axis) and 7t/(l + ft) 
(i/-axis). In the figure the flow lines start on the :r-axis with 7 j = 0. For all of them 7 t = oo at 
some finite value s c , which depends on the initial values. 



In addition, due to the divergence of "ft also Z goes to the strong coupling regime, 
leading to an enhancement of the specific heat, which is however hardly observable in 
two dimensions. Close to the value s c , one has from eq.JIOJ that 

(6.6) 7t «( s - Sc )-\ 
which together with the ea. (|6.3|) for Z gives 

(6.7) Z^(s-s c )~ 3 , Z t « (s-s c r 4 . 
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To convert the above behavior as function of temperature, one may reason in the fol- 
lowing way. In general, the divergence of the length scale corresponds to a vanishing 
temperature, as required by the diffusion law condition L 2 — D/T. In the present case, 
however, the parameter Z renormalizes the temperature so that one has the renormalized 
condition L 2 = D/(ZT). At finite length scale, the vanishing of T is compensated by 
the divergence of Z in such a way that T w (s — s c ) 3 . This implies for the specific heat 
and spin susceptibility 

(6.8) ^r-\ X ~T- 4 /3. 

We may finally notice that the inclusion of the cooperon contribution would modify the 
above group equations, without qualitative changes in the overall behavior (see ref . |74| 
for details ). 

In d = 3 one has a richer scenario depending on the initial value of the running 
variables. In the limit of large 74 and small t the product i 7t obeys the equation 

(6.9) ^ = 

as as as 2 

which has a fixed point for t c "ft,c = £■ This condition gives the asymptotic expression for 
a critical line in the t — j t plane. Close to this critical line, for large values of 7 t, one has 
the approximate solution 

(6.10) t(s) = t(s )e-< s - s «>/ 2 , 

(6.11) 7t " 1 ( S ) = 7t _1 ( s o) - ^ + «e-C-°)/3. 

e e 

One immediately sees from the above equations that for low disorder, t(s ) < e/ 7 t(so), t 
scales to zero and j t scales to a finite value. In the high-disorder regime, t(s ) > e/ 7t (s ), 
7 t diverges at a finite value of s as in the two-dimensional case and t stays finite. In 
fig |21l we report the flow obtained by numerically integrating the RG equations. We 
note, however, that the strong-coupling runaway flow requires to go beyond the one- 
loop approximation we have presented here leaving open the issue whether this proposed 
scenario is realized or not. An approximate treatment of the two- loop correction is 
possible, but its discussion is well outside the scope of this paper. We refer the reader to 
ref. HI. 

Along the critical line or nearby in the low-disorder regime, by converting to a length 
scale via s = — lnTY = — \xi(tD / L 2 ) ~ 21n(L/i), one has that, while t vanishes, the 
conductivity stays finite 

(6.12) a(L)=t- 1 (L)L-\ 



Equation l|6.12ll may be interpreted in terms of a modification of the scaling law (|2.28() 
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0.2 0.4 0.6 0.8 1 

t/(1+t) 

Fig. 21. - The RG flow, in d = 3, in terms of the variable t/(l + t) and 7t/(l +7t). In the figure 
the flow lines start on the x-axis with 7 t = 0. On the rc-axis, there is value t/(l +t) ~ 0.5 below 
which the RG flows to a state with zero t and finite value of 74. The critical line originating 
from this value of t is shown by a thicker line. At large value of 7t this critical line is well 
described by the approximate formula given in the text. For initial larger values of t, the RG 
flow is qualitatively similar to the two-dimensional case. 



due to the scaling behavior of t c ~ L xt , i.e., 

(6.13) n=v{e-x t ). 

On the critical line Xt = e and fj, = 0. On the other hand 74 compensates the vanishing 
of i and diverges like 

(6.14) 7t ~ L £ . 
The equation for Z , in addition, gives now 

(6.15) ^-ijeZ, =^~L 3£ . 

ds 2 



In contrast to the finite value of a, according to the eqs. (|5.45p . ()5.46l l the spin and heat 
diffusion constants vanish, due to the divergenece of Z t and Z . Higher order correction 
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terms could modify this result leading to a vanishing conductivity also in agreement 
with the experiments. Alternatively we can hypothize that the critical line does not 
characterize the metal-insulator transition. The strong spin fluctuations which are the 
relevant physical effect associated with D s — > (x — ► oo) lead instead to an instability 
line before the localization takes place. In this case the system before reaching the 
instability, should make a crossover to one of the universality classes with magnetic 
couplings, which will be discussed later. By assigning scaling exponents in terms of 
inverse length scale xz = — 3e and Xz 2 = ~~ 4e to Z and Z2 and considering that the 
combination ZT ~ l-( x t+xz) mus t scale as Dq 2 ~ L~ 2 since D remains finite, one gets 



At e = 1, one has the temperature power laws —3/5 and —4/5. Hence, as for the two- 
dimensional clear prediction of the theory is a low temperature enhancement of 
the specific heat and spin susceptibility, the latter being generally stronger. A stronger 
enhancement of the spin susceptibility has been indeed observed in Si : P 26, . Theoret- 
ically, as mentioned before, the divergence of the spin susceptibility at low temperature, 
as predicted by the renormalization group flow, has led to the suggestion that the system, 
because of the slowing-down of spin diffusion, tends to form regions of localized magnetic 
moments, which would eventually drive the system into the universality class of mag- 
netic impurities |721 175] . In a number of experimental papers|SU| |27| El OH- where the 
enhancements of specific heat and spin susceptibility have been compared systematically, 
some sort of an effective two-component system made of localized and itinerant electrons 
has been proposed to interpret the data. 

In the presence of any mechanism that inhibits the spin fluctuation enhancement, the 
localizing term in the eq.JOJ f° r t dominates and one has a metal-insulator transition 
with t c ~ 0(e). For instance, in the presence of a magnetic field only the ladders in the 
triplet channel with projection M — ±1 are suppressed. From ea. H5.50fl for the correction 
to the thermodynamic potential, one sees immediately that the spin susceptibility is 
no longer singular, i.e., Z t is invariant upon renormalization. By using ea. H5.57|) and 
eliminating the contribution of the triplet component with M = ±1 in ea. (|6.3|) for Z one 
gets 



(6.16) 



xt = 2 — xz = 2 + 3e 



which yields for the specific heat and spin susceptibility 



(6.17) 



y-3e/(2+3£) ^ y-4e/(2+3e) 



(6.18) 



dZ 
d7 



(6.19) 



ds ~ 2 



The above equations have a fixed point 7 ( * = 1 with a constant Z. It is now direct to 
obtain the equation for the parameter t. After suppressing the M — ±1 triplet channel 
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contributions in ca. H6.lfl and by using the fixed-point condition 7 t * = 1, one obtains 

(6.20) 4^ = ~4 + ( 2 - 2 ln 2)t 2 , 

ds 2 

which has a fixed point t* = e/(2(2 — 2 In 2)) and gives /x = 1 as in the non-magnetic 
impurity case for the non-interacting system. 

In the case of magnetic impurities or spin-orbit scattering, we have seen in subsection 
14 51 that only the ladder in the singlet channel remains diffusive (cf. eas. <|4.38[) and 

Since all triplet channels are massive, 74 drops out in the equations for Z and t, which 
now read 

dt t n 

(6.21) £=- e 2 + * 

dZ t 

^ -dl = -2 Z - 

The equation for t gives the fixed point t* — e/2 and conductivity scaling exponent fi = 1. 
By using the fixed point value for t in the equation for Z ', one obtains 

(6.23) Z«e- (e/4)lns =T £ / 4 . 

Notice that the identical behavior for the cases of magnetic impurities and spin-orbit 
scattering only holds when neglecting the contribution of the pure localization effects. 
When the latter is also taken into account nothing happens for the magnetic impurities 
case, since all cooperon ladders are massive (cf. eas. H4.4Uf) and (|4.41l) '). For the spin-orbit 
case, on the other hand, the ladder in the cooperon singlet channel is still diffusive and 
contributes by minus one half to the standard localization term, as we have discussed 
at the end of subsection 14 51 Hence the combination of the antilocalizing contribution 
from pure interference with the localizing term due to interaction in the singlet channel 
does not change the qualitative behavior of t and gives fj, = 1, even though the fixed 
point value and the approach to it will differ giving Z ~ T £ / 2 . This is relevant in 
the experiments|4"71 l4"51 14^1 151]] (already discussed in points 2), 3), and 4) of subsect. 
I4'6f) . where a value of /i = 1 is observed both in the absence and presence of magnetic 
field. This is exactly what is predicted by the present theory of combined disorder 
and interaction effects, where a magnetic field simply controls the contribution of the 
antilocalizing pure interference effect in the Cooper channel and changes the approach 
to the fixed point. Such a change is indeed observed in From a theoretical point of 
view, we finally comment that in order to perform a quantitative analysis in the presence 
of both diffuson and cooperon diffusive channels requires however the inclusion of the 
interaction amplitude V3 in the Cooper channel. For details we refer the reader to |78l l6]. 
We hence see that all the universality classes share the same conductivity exponent /i = 1 , 
but differ as far as the behavior of Z (and hence of the specific heat) is concerned. To the 
best of our knowledge there are no experiments available to check this last prediction. 
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We are at the end of our journey through the fascinating world of disordered electron 
systems. It is time to draw a few conclusions. We have seen that the non-interacting 
theory is not sufficient to interpret the existing experiments. In two dimensions the most 
relevant result is the prediction of the metallic phase, which is observed to be suppressed 
by the magnetic field. However, a full account of the experimental situation is far from 
being reached. In three dimensions, the predictions of the theory of disordered interacting 
electron systems agree with the experiments whenever there is a magnetic coupling in 
the system and most of the puzzles met while discussing the non-interacting case are 
resolved. In the general case, with no magnetic coupling present, although the strong 
enhancement for the specific heat and spin susceptibility predicted by the theory appear 
to be confirmed by the experiments, a deeper understanding is clearly needed and further 
theoretical and experimental work is required with particular emphasis on the magnetic 
instability problem. 



R. Raimondi acknowledges partial financial support from MIUR under grant COFIN 
2002022534. C. Di Castro acknowledges partial financial support from MIUR under 
grant COFIN 2001023848. 

Appendix A. 
Useful integrals 

In the evaluation of diagrams, as we have seen, we leave the q-integration at the end, 
i.c, the integration over the momenta flowing in the ladder. Due to the presence of the 
diffusive pole, which makes the small-q region dominant, the remaining integrals over 
the fast momenta can be expanded in powers of q and ui. To this end it is useful to 
expand the Green's function as 



where on the right-hand side G = G(p, e). Then, the integration over the fast momenta 
involves integrals containing products of retarded and advanced Green's functions with 
the same argument. Finally, by using the residue theorem and the formula for the residue 
of poles we get the useful formula 



* * * 



(A.l) 



G(p + q, e + u) = G [l - (w - v • q)G + (w - v • q) 2 G 2 + ...] , 



(A.2) (m,n) = J2(G R r(G A ) 



P 



n 



(-1) 



m i n+m 2irN< 



o 



n(n + l)...(n + to — 1) 
(m-1)! 



.n+m— 1 



The most frequent cases are 



(A.3) 
(A.4) 
(A.5) 



(1. 1) = 27rJV T 
(2, 2) = 2 2irN T 3 

(1.2) = i 2ttA t 2 
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(A.6) 


(2, 1) = -i 2ttN t 2 


(A.7) 


(1,3) - - 2ttN q t 3 


(A.8) 


(3, 1) = - 2^V r 3 


(A.9) 


(1,4) = -i 2ttN t 4 


(A.10) 


(4, 1) = i 27tA t 4 


(A.ll) 


(2,4) = -4 2ttN t 5 


(A.12) 


(4, 2) = -4 27dV T 5 


(A.13) 


(3,3) = 6 2itN t 5 . 



Integrals containing scalar products are evaluated as 

2 2 

(A.14) ^(p ■ <l) 2 (G R ) m (G A ) n = ^j- J2(G R ) m {G A ) r 



Appendix B. 

Magnetic impurities 

To see this, let us consider the following term in the Hamiltonian 



H disorder = d r ipl{r) [u(r)S a/ 3 + u s (r)S ■ er Q/3 ] tpp(r), 



(B.l) 



where S is the spin of a paramagnetic impurity located at r having a scattering amplitude 
U s (r), whereas u(r) is the scattering amplitude due to non magnetic impurities and 
already taken into account within the self-consistent Born approximation. The first step 
amounts to recompute the Green's function in the presence of the full term ea. (|B.ljl . 
One gets for the Green's function 



(B.2) 



G H (p,e) 



where r 1 = 27rA r u 2 (r), r s 1 = 2ttNqu 2 s (t)S(S + 1) and for the single impurity line 

1 



(B.3) 



Uct/3j5 



2ttN q t 



Suf3^5 + O a ft ■ CT 7< 5 

3r s 



where the meaning of the spin greek indices is shown in fig l22l To find the expression for 
the diffuson and cooperon ladders in the presence of magnetic impurities, it is convenient 
to exploit the conservation of the total spin of the pair (particle-hole for the diffuson 
(a — S) and particle-particle for the cooperon (a + 7)). To this end we use the singlet 
and triplet spin projection operator, which for the diffuson are given by 



(B.4) 



1 
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a 



Pi 



P 4 



Pi 



a 



P 



P 2 



P 3 



P 



P 2 



"PS y > ^ 5 "P 4 

Fig. 22. - A single impurity line with the spin and momentum structure which are related to 
the diffuson and cooperon ladder. (See fig |121 . 



(B.5) Tapys — 2 <JaS ' a Pi = 4 fifiapSys - • a^s] , 



and for the cooperon by inverting one of the Green's function: 



( B -7) T p ~ p s = - [3S a ^ s + o a p ■ cr~<5} 



where S and T stand for singlet and triplet, respectively. If we indicate with the L^ S,T 
and Lc the the single-impurity line contribution for both singlet (S) and triplet (T) 
components, 

1 f i + L 



(B.8) 


L (0)S 


= Tr{US p ~ h ) 


(B.9) 


L (0)T 


= Tr(UT p - h ) 


(B.10) 


L (0) S 


= Tr(US p ~ p ) 


(B.ll) 


L w 


= Tr(UT p ~ p ) 



2ttN t \ T, 

1 'l ' 



2ttN t \ 3r 

1 U r 



2ttN t \ T, 

1 * 



2ttN t V 3r 
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one has for L S ' T and Lf ' T the diffuson and cooperon ladders 

r(0)S,T r(0)S,T 



i-nL(°) s . T ' c l-n4 0)S ' T ' 

The trace sign implies summation over four spin indices. The quantity II is defined as 
in ea. l4.11|) with the appropriate Green's function for this case. Its evalution leads to 
ea. pTTS)) with 1/r replaced by 1/t + l/r s sw 1/t(1 - r/r s ). 

Appendix C. 

Spin-orbit scattering 

In the presence of spin-orbit interaction, in the scattering Hamiltonian one adds a 
term 

(CI) H so = — ^ E ^(r) W a p ■ W(r - R) A p] ^(r) 

R 

where R indicates a ion site and V(r — R) is the corresponding potential. The matrix 
element between states of momentum p and p' is 

(C.2) < p\H so \p' > = -1-2 E a ^ ■ [ dr e ~ iP "~ t W ( r - R ) A P] eiP ' r 

c j 



R 



(C.3) = iu 



R P ^ 



Since the scattering depends on the momenta of the particles involved, one has to 
consider separately the contribution to the self-energy and to the single-impurity line in 
a ladder resummation. For the self-energy one has 

m a\ ^ i s ~~o (p A p') (p'Ap) . . . 

(C.4) E SOQ/3 (p, e) = -u z so 2 CT «7 2 ' CT 7/3 G (P > e ) 

p , Pf Pf 

fnn\ .sign{e) 
(C.5) = -i— 5 af 3 

so that the Green's function reads 



(C6) G R (p,e) 



-(- + — 

2 \T T, 



For the contribution of a single impurity line one has, by performing the angle average 
(indicated with a bar) over p and p', 



rr"7\ tt i \ ~2 (Pi A P2) (P3Ap 4 ) 

(\>- I ) Ua/3~/6{Pl, P2, P3,Pi) = — U so 2 ' CJa l 3 2 ' °l f > 

Pf Pf 
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with the condition that the total momentum is conserved, Pi +P3 = P2 + P4- In the case 
of the diffuson ladder one has pi ~ P4, while for the cooperon p x ~ — p 3 (See flg l22)l . 
As a result 

(C8) ^ = _1_^.^ 

(C9) uP-P s = --l— aa ,. ajS . 

As in the case of magnetic impurities one calculates the single-impurity line contribution 
(CIO) L^ s = -1— + — 

(di) i (0)T = ^r-(i T 



2ttN t \ dr. 

(C13) 4° )T = -i— ( 1 7 



2ttN t V dr. 

In the diagram giving the weak localization correction, in general the spin structure has 
the form 6 a s6yp. By using the projection operators developed in the previous Appendix, 
one can separate the singlet and triplet contribution to the weak-localization correction 
in the Cooper channel 

(C14) Tr(S p a - p lS S aS S 7l3 ) = -1 

(C15) Tr(T^ s S a5 S j!3 ) = 3. 

We see that, while the triplet is localizing, the singlet has an opposite effect. In the 
absence of spin-orbit scattering both the singlet and triplet are massless and since the 
triplet contribution is three times larger, its effect prevails. In the presence of spin- 
orbit scattering, the triplet becomes massive and does not contribute to the logarithmic 
singularity term in d = 2. The latter therefore changes sign and becomes antilocalizing. 



Appendix D. 

The long-range interaction case 

The effective screened interaction is given by 

V c (q) 2ne 2 Dq 2 -iuj 



(D.f) V H (q,u;) 



l + V c {q)K ot) (q,uj) q Dq 2 + Dqn - iuj ' 



where we used the two-dimensional expressions for the Vc = 2ire 2 jq and the Thomas- 
Fermi inverse screening length k = 47re 2 Ao. In the energy and momentum region given 
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Fig. 23. - Diagrams for the conductivity to the lowest order in the interaction, (a) and (b) 
Exchange, (c) and (d) Hartree. Interaction is shown as a thick dashed line. 



by u) < Dqn and Dq 2 < u>, the correction to the density of states becomes 

1 f°° 1 
(D.2) 5N(e) = —lm / doj /(« - e) £ 

bj/DK<\q\<yfuj/D 

(D.3) = -* M ,e| T )l„(J--). 

We see that the correction becomes log-square! This is a peculiar feature of the single- 
particle density of states. In fact, all other physical quantities that we deal with in 
these lecture notes acquire logarithmic corrections only even in the presence of Coulomb 
interaction. As we shall see in Appendix iHl the density of states can be reabsorbed into 
the definition of the scattering amplitudes and drops out from the renormalization group 
equations. 

Appendix E. 

Details on the evaluation of the interaction correction to the conductivity 

In this Appendix, we show how to obtain the correction to the electrical conductivity 
due to the combined effect of disorder and interaction. The diagrams contributing to 
electrical conductivity to lowest order in the interaction are shown in fig l23l Diagrams 
(a) and (d) are obtained by inserting a self-energy correction into the Green's function 
and one has to consider also the symmetric ones with the self-energy insertion in the 
bottom electron line. Diagrams (b) and (c) are due to vertex corrections. We begin 
our discussion with diagrams (a) and (b). The extension to diagrams (c) and (d) is 
straightforward. The expression for diagram (a) reads 

^ Q) (r,r';r») = 2T^ 7 l (r) 7 ^(r') f dr x dr 2 T^ V(n, r 2 ;w m ) 
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(E.l) 



<5(r, ri; e„ + fi„)£?(ri + Q u - w m )£(r 2 , r'; e. 



+ 0„)a(r',r; e, 



)• 



To this expression one has to add the one corresponding to having the interaction line in 
the bottom electron line. 7 l (r) is the real-space representation of the current vertices. We 
do not need to write here its explicit expression since at the end, after the impurity av- 
erage, we recover translational invariance and go back to the momentum representation. 
Diagram (b) is instead given by 



iig ) (r,r / ;0)=2T^7 i (r) 7 i(r') / dridr^r, r x ; e n + n v )Q{v 2 , r; e n ) 



In both the eas. (|E. 1( 1 . (|E.2(1 . the Matsubara frequency sums are transformed to contour 
integrals in the complex plane by means of standard manipulations. Then one gets for 
diagram (a) (including the other diagram with the two electron lines interchanged) 



f(e - n)b(u)(V R - V A )^G?G?^G?{Gtn ~ G A _n) 
+ f(e - - e)V R G R (G R _^ ~ Gf_ JGf (G?_ n - G?_n) 

+ f(e)b(u,)(V R - V A ) u G R G R _ u G?Gtn 
+ /(e)/( w - e)V R G R (G R ^ - G A _JG R G?_ n 

- f{e)b{u){V R - V A )^G A Gt.G A Gtn 

- /( e )/( w - e)V A G A (G R ^ - GtJGfGtn 

+ f(e Q)b(")(V R - V A ) u G R G R _ u _ n Gf_ n G R _ n 

+ f(e - - e + W£G R G R _^Gt^ ~ G A _ w _ n )G R _ n 

- f(e - n)6( w )(V* - V A )^GfGtnGt^Gtn 

- f(e - il)f(u - e + n)^ A Gf Gt_ n (G*_ u _ n - G A _JG A _n 
+ f(e)b(u)(V R - V A UG R - G A )GtnGt.-nGtn 

(E.3) + f(e)f(u, - e + n)V A (G R - Gt)G£_ a (G?_ u _ a - Gt„-n)Gtn] • 



(E.2) 



T£0(: 





Diagram (b) gives 




GtJ 
)-n-u)G e - 
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- f(e)f(u - e - W<?G*. n Gi(G*. n _ u - G? +n _ u )G?_ u 
+ f(e + Q)b{u)(V R - V^G^GtG^Gt. 

+ f(e + - e)V*G? +n G?G? +n ^(G^ ~ C£_J 

+ f(e + - e - fi)V« G? +n G?(G? +n _„ - G? +n _JGtu 

- f(e + - V A ) u Gt^GfGt +a _ u Gt_ u 

(E.4) - /(e + - e - 0)V»Gi +a Gt(G*. n _ u - G? +n _jGt u _ 



Since we are discussing the frequency structure, in eas. (|E.3|) and (|E.4|) we have 
dropped the explicit dependence on space coordinates. The Green's functions are pre- 
sented in the same order as in eas. (|E.l[i and l|E.2(l where the space dependence is shown. 
In eas. (|E.3(l and (|E.4|) we may now perform the impurity average. First, averaging im- 
purity pairs belonging to the same Green's function line implies the replacement of the 
Green's function with its self-consistent Born approximation expression. Secondly, we 
have to perform the average of impurity pairs belonging to different Green's function 
lines. This can be performed by arranging the Green's functions on the sides of a square 
and inserting ladders wherever possible. At the leading order in the expansion parame- 
ter, we neglect, as a rule, all the diagrams in which a crossing of impurity lines occurs. 
Depending on the sequence of retarded and advanced Green's functions around the sides 
of the square, we may insert two or three ladders. For instance, terms with four retarded 
Green's functions give zero since all poles lie on the same side of the real axis. In ea. (|E.3|l . 
the terms that allow two or three ladder insertions are the second, fourth, sixth, cigthth, 
tenth, and twelfth: 

+ f(e - Cl)f( u - e)V*(q)j^G?(G?_ u ~ Gt_JG*(G*. a - G A _ Q ) 



+ /(e)/(w - e)Kf(q)7* 7 iG«(G!f_ w - Gf_ u )G* G e _ n 



/( e )/(w - e)V u A (q)y 7 iGi(G?_ u - Gt_ u )G?G£_ 



+ f(e - fi)/(w - e + n)V u R (d)r 7 J G«Gf_ n (Gf_ w _ n - Gt^n)Gtn 



f(e - fi)/( w - e + n)V^G*Gj_ a (G R _ u _ n Gf_JG^_ n 
(E.5) + /(e)/(w - e + (q)r 1*{G R ~ G A )G A _ n (G R _ n - G?_ u _ n )G A _ 



In the same way, from ea. HE.4f> we pick up the third and the eleventh term 



+ /(e)/(w - e - W?^)l l l ] G R G?\G R u - G A JG A _ 



(E.6) - /(e + Q)f(u - e)V A (q)Yj>G A G A G R JG^ - G A _J 
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(b) 

to 



Fig. 24. - Effective conductivity diagrams. 



The b ar ove r the prod ucts of Green's function indicates the impurity average. Notice that 
in eas. (|E.5J) and (|E.6J) the current vertices appear under the impurity average bar. After 
the average and the restoration of traslational invariance, there appear summations over 
momenta. The sum over the slow momenta that enter the interaction and the ladders 
are performed at the end, while the sum over the fast momenta entering the Green's 
functions are performed with the help of residue theorem within the approximations 
explained in detail in Appendix To this end, all frequencies in the Green's function 
can be set to zero in the leading order in the diffusive regime expansion lot < 1. As a 
result, the impurity average of the product of Green's functions does not depend on the 
energy. Since, the ladders depend only on the slow frequencies uj and fl, we can perform 
the e-integration at once by using the useful identity 

(E.7) |~ d e /(e) f{u - e) = | (coth (^) - l) = F(u). 

We may finally consider explicitly the impurity average. To illustrate the procedure, 
let us consider the first term in ea. (|E.5ll . It contains four products of four Green's 
functions each. The first product GfGf_ M GfGf_ q gives zero upon averaging. The term 
G^G^^Gf Gf_ n vanishes because of the vector nature of the current vertices. There 
remain the terms GfG^_ u G^G^_ n and GfGf_ u Gf-Gf_ n . Upon averaging, the first 
term gives rise to an effective diagram with three ladders, corresponding to (a) of figl2*4l 
The second term, on the other hand, yields effective diagrams with two ladders only, 
corresponding to (c) and (d) in fig l24l By following this line of reasoning, one obtains 
all the diagrams of figE] including those obtained by interchanging the top and bottom 
Green's function lines. One key point to notice is that the diagrams with two ladders 
cancel each other, i.e., the sum of (c),(d), and (e). This cancellation is shown in detail in 
rcf. 59 . To sec it, let us consider the terms in eas. ljE.3IE.4ll which contain the retarded 
interaction V^f (a similar analysis can be done for the terms containing V^y). By using 
ea. (|E.7jl . we get from eas. (|E.5IE.6|) 
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-Fiuj^GRGt^Gtn ~ F^h^GfGt^G^ 
-F{" n) 1 * 1 iG? +n G?G£, n _ a GtJ) , 

which is readily seen to vanish by considering the following averages: 



l^G^Gt^Gf_ 



l^G^Gt^Gt 



^G^ +a G^Gf +n _ u Gt, 

As a final comment, we note that each individual diagram with two ladders, in the 
presence of long-range forces, will suffer from the strong singularity as in the case of 
the density of states, as discussed in Appendix \B\ However, it has been shown that 
the singularity due to the long-range forces can be incorporated into a gauge factor, 
which drops out in the evaluation of gauge-invariant quantities. This is the origin of the 
cancellation of the diagrams with two ladders [75], 

Let us now turn our attention to the diagrams with three ladders. The impurity 
average needed for the diagram(a) of figl2His 

^ 7 ! Gf(p)e.(p)G^(p-q) 

p 

$> J 'Gf (p')GtMGtn(p' q)(r°(q,c)) 2 L(q,c- O) 
p' 

(-4TreN DT 2 q i ){-4'KeNoDT 2 q j ) 1 

2ttN t 4 (Dq 2 -ia;) 2 (£>q 2 -i(a;-fi)) 

. £W 

77(70 (L>q 2 - kj) 2 (L>q 2 - i(w - SI)) 



(E.8) 



where the factors in round brackets, in the second line, arise from the integration of 
the three Green's functions with a current vertex. Notice that in diagram (b) of fig l24l 
the two integrations over products of three Green's functions produce an opposite sign. 
This gives an overall minus sign for diagram (b) with respect to (a). In the last line 
(T = 2e 2 N D is the Drude conductivity. By collecting in eas. i|E.5|l (first, third, fourth, 
and sixth line) and l|E.6(l all the terms giving rise to diagrams with three ladders we get 



R i i(0,Sl) = -4a J2 I 
q J ~ 

F{u - O)Kf(q) 



2^ 



Dq*q> 



(Dq 2 - iw) 2 (L>q 2 - i(w - ft)) 



(Z)q 2 + iw) 2 (L>q 2 + i(w - SI)) 
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{Dq 2 -iuj) 2 (Dq 2 - i(w + ft)) 



(E.9) 



+ F(u; + n)V A (q) 

- f(lu - n)C(q) 

- F(w - n)C(q) 

- F(u> + fl)V A (q) 

- F(w + n)Kf(q) 





(Dq 2 


+ iw) 2 (Dq 2 + 
Dq l qi 


i(u> 


+ 0)) 




~\oj) 2 {Dq 2 - 
Dq l qi 


i(u> 


+ «)) 




-iuj) 2 {Dq 2 - 
DgV 


i(u> 


-0)) 


(Dq 2 


+ \uj) 2 (Dq 2 + 
Dq l qi 


i(u> 


+ n)) 


(Dq 2 


+ \u) 2 (Dq 2 + 


i(u> 


-0)) 



Notice that the first term cancels with the sixth and the fourth with the seventh. By 
diving by — iQ and sending fl to zero we get ea. (|5.11|l quoted in the text. 

Appendix F. 

Details of the evaluation of the thermodynamic potential 

To begin with, let us consider, at fixed impurity configuration, the first-order exchange 
interaction correction to the thermodynamic potential (see the second diagram in fig. 1161 
without the inserted ladder) 



(F.l) Afi = j ( 1 y^ 2 E ydrdr'V(r-r', Wm ,A)^(r,r'; e „)g(r',r;e ri 



where we have used the standard trick |44| to multiply the interaction by a parameter 
< A < 1 

V(r - r', u m , A) = AV(r - r', uj m ). 
The sum over the Fermionic Matsubara frequency gives 

T^2g(r,r';e„)g(r',r;e n +Lj m )\ iuJm=UJ+i0+ = 

£ n 

' oc de 



2ni 



/(e) [(G>, r'; e) - G (r, r'; e))G R (r', r; e + w) 



(F.2) 



+ G A (r, r'; e - u>)(G R (r', r; e) - G A (r', r; e))] 



By performing the impurity average we need to keep only terms with both retarded and 
advanced Green's functions. This selects in ea. l|F.2(l the term 
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= r £ [/(e +uj) ~ /(e)] GA(r ' r,; e)Gfl(r '' r; e + u) 

(F.3) - ^(2^ T) 2 L(q, W ), 

where, making use of the fact that the impurity average of the product of Green's func- 
tions does not depend on the energy, we have performed the integration over the fre- 
quency. As a result, ea. (|F.l|> becomes 

1 ^t'W N V(q,u} m , \)\u> m \ 
A Dtf + l^l ' 

from which, by taking iVoV(q, ui m , A) = X(Vi — 2V2), one obtains the ea. l|5.14|l of the 
main text. Whence the dynamical amplitude resummation is inserted into the diagrams 
of fig^jjone obtains the final expressions l|5. 51)11 . 1)5. 511) for the thermodynamic potential. 



(F.4) 



An = - 



Appendix G. 



Ladder in the presence of Zeeman coupling 

In the presence of a magnetic field, electron energies are changed by the Zeeman 
energy, so that the Green's function reads 



(G.l) 



G£(p,e) = 



1 

27 



where uj s — q^bB with g the Lande factor and fis the Bohr magneton. The spin 
projection takes values a = ±1/2. One sees that the Zeeman energy au s enters the 
Green's function as an energy in shift. This allows to get immediately the ladder in 
the presence of magnetic Zeeman coupling, since the energy difference to is shifted by 
the difference of the Zeeman energies of the particle-hole pair. For instance, by making 
reference to the spin structure of fig l22l and taking into account the spin conservation 
along a Green's function line, one has 



(G.2) 



L aS {q,aj) 



1 



1 



2ttN t 2 Dq 2 - iuj - i(a - S)uj s 



This show that only the triplet components with total spin projection a — S = M = ±1 
are affected by the magnetic field. 



Appendix H. 



The ladder renormalization 



In this Appendix we show that the logarithmic corrections found for the physical 
quantities can be absorbed into a renormalization of the parameters characterizing the 
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e+co-£2 
P 

(a) 



V V 



(b) 




(c) 




Fig. 25. - Diagrams for the ladder self-energy. A similar set of diagrams is generated by the 
interchange of the interaction line (tick-dashed line) between top and bottom Green's function 
lines. 



ladder propagator. This identification is the formal basis of the renormalizability of 
the effective field theory, whose physical meaning is discussed in the text via the Ward 
identities. Let us assume that the ladder, in the presence of interaction, gets renormalized 
as 



(H.l) 



1 



1 



^ 2 



2ttN t 2 Dq 2 



2ttN t 2 D R q 2 - iZuo 



where £, D Rl and Z represent the effective wave function renormalization, the renormal- 
ized diffusion coefficient, and the renormalization of the frequency. By expanding 



(h.2; 



i 



(2SC - SD/D)Dq 2 - i(25( - SZ)lu _ T, L {q, u>) 



D R q 2 - \Zuj Dq 2 



(Dq 2 



{Dq 2 



and the last equation defines the ladder self-energy. The diagrams for the self-energy are 
shown in fig l25l The first step is the integration over the fast momenta running within 
the Green's functions. This integration amounts to the evaluation of several integrals of 
products of the type (G R ) m (G A ) n , whose result is given in Appendix lAl For diagrams 
(a-c), in the small k, O, q and u> limit, we obtain 



(H.3) 



I abc = (2ttN t) 2 2ttN q t 4 [D(q 2 + k 2 ) - \{w + fi) - 2q • k] , 



where the first factor (2ttNqt) 2 represents the two integrations over the two Green's 
functions at the interaction vertices. The rest gives the integration over the remaining 
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Green's functions. In a similar way, integration over Green's functions gives for diagram 
(d) 

(H.4) I d = -(2itN q t 2 ) 3 . 

The diagrams (a),(b),(c) shown in fie 1251 finally yield 

(k, q, ft m , w„)V(k, f2 m ) 

( 5) ^,a6c--22 2. (2^V T2)3pp + |a„|) 2 

and diagram (d) 

i d v(k, n m ) 



(H.6) SL, d = 2T J2 E 



, , k (2^iV T 2 )3(^(k + q) 2 + \n m + <j v \) ' 

In eas. ljH.5IH.6|l we have used Matsubara frequencies. The relative minus sign comes 
from the integration over the fast momenta. The factor of 2 is due to the fact that there 
is another set of diagrams generated by interchanging the interaction line between the 
two electron lines. One may check that the sum of eas. ljH.5IH.6p vanishes in the limit 
q = and u) = 0. For small, but finite external momentum and frequency, we rewrite 
cq. (|H.5|l in the form 

^L,abc --21 ^ 2^ Dk 2 + |n| 

f\ 0T ^ (Dg 2 + \^\)v(k,n m ) 



and ea. (|H.6|l as 

^ = 2T £ £ v(k '" m) 



, £>(k + q) a + |fi m + o;„| 

e ri <il r(l <e n -\-LO l/ k 



(H.8) +2T £ £ 



V(k,Q ro ) V(k,Q m ) V(k,Q r , 



Dk 2 + |ft m | D(k + q) 2 + |r! m +^| f?k 2 + |a 



where less divergent terms have been neglected. The first term in the square brakets 
of ea. i|H.8|) cancels with the first term in ea. (|H.7|) . Let us analyze the first term of 
eq. (|H.8fl . By transforming the Matsubara sum into an integral in the complex plane and 
analytically continuing e — > —ie and e + uj — > — i(e + u), it reads 



T 1 = IT V V V ( k > n rn) 

L - 4 so,+ ^ D(k + c0 2 + \n m +u; v \ 



2 f 00 

-— J dfi[/(n-e-w)-/(n- e )]X; 



F fl (0,0) 



L>k 2 - + w) 
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2lo x 



V R (0,0) 



2m V -Ck 2 - i( e + w ) 

k 



. - 2V 2 . / 1 

la; rr— — In — 

4tt 2 L> V er 



(H.9) 



-iw I 3 



having used f(Cl — e — u>) — f(Cl — e) sw uj — e). In the last line we have also included 

the contribution of the Hartree diagrams. The difference between the second and third 
term in the square brakets of ea. ljH.8[l may be expanded in powers of q and ui v . The 
lowest order term reads 



2T E E v ( k '^< 



2D 2 q 2 \i 2 



{Dq 2 + M) 



(Dk 2 + \n m \r (Dk 2 + \n m \) 2 

Dq 2 Dk 2 V R (0, 0) (Dq 2 - ioj)V R (Q, 0) 



(Dk 2 - ifi) 3 



(L>k 2 - ifl) 2 



y 4tt 2 L> V er 



(H.10) = L»g 2 I 2 - (iV - iw)Ii. 



fn 2 • /i-2^2 , 
(L><7 — iuj) — . „ — In 



4tt 2 L> 



Finally, the second term in ea. l)H.7(l reads 



^L,abc = -(Dq 2 + \lu v \)2T J2 E 



v(k,n r . 



(dp + |n m |) 2 



-(A? 2 - ho 



2ni 



dO/(f2 



y fl (o,o) 

(D/c 2 - ift) 2 



47T 2 L» 



In I - 

er 



(H.ll) = -(Dg 2 - ho) Ii. 

One sees that the by inserting the self-energy results ljH.9IH.ll|) into ca. (|H.2|l one gets 



(H.12) 
(H.13) 
(H.14) 



C = i -ii, 

Dr 
D 

Z = l- I 3 . 



= 1-Ia 



In the main text the renormalized diffusion coefficient Dr will be renamed D. We con- 
clude this Appendix by remarking that the above renormalizations coincide with the 
perturbative corrections of the single-particle density of states, ea. l|5.1U|l . conductivity, 
ea. (|5.12|l and specific heat, ea. l|5.17|) . satisfying at this order the Ward identities identi- 
fications. 



Disordered Electron Systems 



71 



REFERENCES 

[1] BERGMANN G., Phys. Rep., 107 (1984) 1. 

[2] Lee P. A. and Ramakrishnan T. V., Rev. Mod. Phys., 57 (1985) 287. 

[3] Altshuler B. L. and Aronov A. G., in Electron- Electron Interactions in Disordered 

Systems, edited by Pollak M. and Efros A. L. (North-Holland, Amsterdam) 1985, 

p. 1. 

[4] Castellani C. and Dl Castro C, in Lecture Notes in Physics 216 edited by 

Garrido L. (Springer- Verlag) 1985. 
[5] Castellani C, Di Castro C. and Strinati G., in Lecture Notes in Physics 268 edited 

by Garrido L. (Springer- Verlag) 1985. 
[6] Finkelstein A. M. , Soviet Scientific Reviews, 14 (1990) . 
[7] Kramer B. and MacKinnon A., Rep. Prog. Phys., 56 (1993) 1469 . 
[8] Belitz D. and KlRKPATRICK F. R., Rev. Mod. Phys., 66 (1994) 261. 
[9] Ioffe A. F. and Regel A. R., Prog. Semicond., 4 (1960) 237. 
[10] Mott N. ¥.,Philos. Magazine, 26 (1972) 1015. 

[11] Dolan G. J. and Osheroff D. D., Phys. Rev. Lett., 43 (1979) 721. 

[12] Rosenbaum T. F., Milligan R. F., Paalanen M. A., Thomas G. A., Bhatt R. N., 

and Lin W., Phys. Rev. B , 27 (1983) 7509. 
[13] Rosenbaum T. F. Andres K., Thomas G. A., and Bhatt R. N., Phys. Rev. Lett, 

45 (1980) 1723. 

[14] Stupp H., Hornung M., Lakner M., Madel O., and v.Lohneysen H., 

Phys. Rev. Lett, 71 (1993) 2634. 
[15] Rosenbaum T. F., Thomas G. A., and Paalanen M. A., , Phys. Rev. Lett., 72 

(1994) 2121. 

[16] Stupp H., Hornung M., Lakner M., Madel O., and v.Lohneysen H., 

Phys. Rev. Lett, 72 (1994) 2122. 
[17] Castner T. G., Phys. Rev. Lett., 73 (1994) 3600. 

[18] Stupp H., Hornung M., Lakner M., Madel O., and v.Lohneysen H., 

Phys. Rev. Lett, 73 (1994) 3601. 
[19] Thomas G. A., Ootuka Y., Katsumoto S., Kobayashi S., and Sasaki S., 

Phys. Rev. B, 25 (1982) 4288. 
[20] Hertel G., Bishop D. J., Spencer E. G., Rowell J. M., and Dynes R. C, 

Phys. Rev. Lett. , 50 (1983) 743. 
[21] Yamaguchi M., Nishida N., Furubayashi T., Morigaki K., Ishimoto H., and 

Ono K., Physica B (Amsterdam), 118 (1983) 694. 
[22] Rhode M. and Micklitz H., Phys. Rev. B, 36 (1987) 7572. 
[23] McMillan W. L. and Mockel J., Phys. Rev. Lett. , 46 (1981) 556. 
[24] Kobayashi S., Ikeata S., Kobayashi S., and Sasaki W., Solid State 

Communications, 32 (1979) 1174. 
[25] Thomas G. A., Ootuka Y., Kobayashi S., and Sasaki W., Phys. Rev. B, 24 

(1981) 4886. 

[26] Paalanen M. A., Graebner J. E., Bhatt R. N., and Sachdev S.,(1988) 

Phys. Rev. Lett, 61 (1988) 597. 
[27] Lakner M. and v.Lohneysen H., Phys. Rev. Lett, 63 (1989) 648. 
[28] Ikeata S. and Kobayashi S., Solid State Communications, 56 (1985) 607. 



72 



Carlo Di Castro Roberto Raimondi 



[29] Paalanen M. A., Sachdev S., Bhatt R. N., and Ruckenstein A. E., Phys. Rev. Lett, 
57 (1986) 2061. 

[30] Alloul H. and Dellouve P., Phys. Rev. Lett., 59 (1987) 578. 

[31] Hirsch M. J. Holcomb D. F., Bhatt R. N., and Paalanen M. A. , Phys. Rev. Lett, 
68 (1992) 1418. 

[32] Schlager H. G. and v. Lohneysen H., Europhys. Lett, 40 (1997) 661. 

[33] Kravchenko S. V., Mason W. E., Bower G. E., Furneaux J. E., Pudalov V. M. , 

and D'lORlO M., Phys. Rev. B, 51 (1995) 7038. 
[34] Abrahams E., Kravchenko S. V., and Sarachik M. P., Rev. Mod. Phys., 73 (2001) 

251. 

[35] Kravchenko S. V., and Sarachik M. P., Rep. Prog. Phys., 67 (2004) 1. 
[36] Anderson P. W., Physical Review , 109 (1958) 1492. 

[37] Mott N. F., Advances in Physics (Philosophical Magazine Supplement, 16 (1967) 49. 

[38] Edwards J. T. and Thouless D. .J, J. Phys. C, 5 (1972) 807. 

[39] Licciardello D. C. and THOULESS D. .J, Phys. Rev. Lett, 35 (1975) 1475. 

[40] Abrahams E., Anderson P. W., LicciardelloD. C., and Ramakrishnan T. V. , 

Phys. Rev. Lett, 42 (1979) 673. 
[41] Wegner F.,Z. Phys. B, 25 (1976) 327. 

[42] Lancer J. S. and Neal T., Phys. Rev. Lett, 16 (1966) 984 . 

[43] Castellani C, Di Castro C, Forgacs G., and Tabet E., J. Phys. C: Solid State 
Phys., 16 (1983) 159. 

[44] A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinski, Methods of Quantum 

Field Theory in Statistical Physics, Dover Publications Inc, New York (1975), Sec 1.3. 
[45] Zala G., NAROZHNY B. N., and Aleiner I. L. , Physical Review B, 64 (2001) 214204. 
[46] Gorkov L. P., Larkin A. I., and Khmelnitskii D. E., Pis'ma Zh. Eksp. Teor. Fix., 

30 (1979) 248 [JETP Lett., 30 (1979) 228]. 
[47] Okuma S., Komori F., and Kobayashi S., in Anderson Localization. Proceedings 

of the International Symposium edited by Ando, T.; Fukuyama, H. (Springer- Verlag) 

1988, p.78. 

[48] Nishida N., Furubayashi T. and YAMAGUCHI M., Solid State Physics, 19 (1984) 
410. 

[49] Katsumoto S., Komori F. Sano N., and Kobayashi S., Journal of the Physical 

Society of Japan, 56 (1987) 2259. 
[50] Katsumoto S., in Anderson Localization. Proceedings of the International Symposium 

edited by Ando, T.; Fukuyama, H. (Springer- Verlag) 1988, p.45. 
[51] Paalanen M. A., Rosenbaum T. F., Thomas G. A., Bhatt R. N., and Lin W., 

Phys. Rev. Lett. , 48 (1982) 1284. 
[52] Shafarman W. N., Koon D. W., and Castner T. G., Phys. Rev. 5,40 (1989) 1216. 
[53] Dai P., Zhang Y., and Sarachik M.P., Phys. Rev. Lett, 67 (1991) 1914. 
[54] Dai P., Zhang Y., Bogdanovich S., and Sarachik M.P., Phys. Rev. B, 48 (1993) 

4941. 

[55] Attaccalite C., Moroni S., Gori-Giorgi P., and Bachelet G.B. , Phys. Rev. Lett, 
88 (2002) 256601. 

[56] Simonian D. Kravchenko S. V., and Sarachik M. , Phys. Rev. Lett., 79 (1997) 
2304. 

[57] Altshuler B. L. and Aronov A. G., Pis'ma Zh. Eksp. Teor. Fiz., 30 (1979) 514 [JETP 
Lett, 30 (1979) 482]. 

[58] Altshuler B. L., Aronov A. G., and Lee P. A. , Phys. Rev. Lett., 44 (1980) 1288. 
[59] Altshuler B. L., Khmelnitskii D. E., Larkin A. L, and Lee P. A. , Phys. Rev. B, 
22 (1980) 5142. 



Disordered Electron Systems 



73 



[60] Altshuler B. L. and Aronov A. G., Solid State Communications, 46 (1983) 429. 
[61] Finkelstein A. M. , Zh. Eksp. Teor. Fiz., 84 (1983) 168 [Sov. Phys. JETP , 57 (1983) 
97]. 

[62] Castellani C, Di Castro C, Lee P. A., and Ma M., Phys. Rev. B, 30 (1984) 527. 
[63] Castellani C. and Di Castro C, Phys. Rev. B, 34 (1986) 5935. 
[64] Castellani C, Kotliar C, and Lee P. A., Phys. Rev. Lett, 59 (1987) 323. 
[65] Castellani C, Di Castro C, Forgacs C, and Tabet E., Nucler Phys. B, 225 
(1983) 441. 

[66] Altshuler B. L., Aronov A. C, and Zyuzin A. Yu. , Zh. Eksp. Teor. Fiz., 84 

(1983) 1537 [Sov. Phys. JETP , 57 (1983) 889]. 
[67] Castellani C, Di Castro C, Lee P. A., Ma M., Sorella S., and Tabet E., 

Phys. Rev. B, 33 (1986) 6169. 
[68] Raimondi R., Castellani C, and Di Castro C, Physical Review B, 42 (1990) 4724. 
[69] Castellani C, Di Castro C, Kotliar C, Lee P. A., and Strinati C, Physical 

Review Letters, 59 (1987) 477. 
[70] Castellani C, Di Castro C, Kotliar C, Lee P. A., and Strinati C, Physical 

Review B, 37 (1988) 9046. 
[71] Di Castro C, in Anderson Localization. Proceedings of the International Symposium 

edited by Ando, T.; Fukuyama, H. (Springer- Verlag) 1988, p.96. 
[72] Finkelstein A. M. , Z. Phys. B, 56 (1984) 189. 

[73] Castellani C, Di Castro C, Lee P. A., Ma M., Sorella S., and Tabet E., 

Phys. Rev. B, 30 (1984) 1596. 
[74] Castellani C, Di Castro C, and Lee P. A., Phys. Rev. B, 57 (1998) 9381. 
[75] Punnoose A. and Finkel'stein A. M., Phys. Rev. Lett, 88 (2002) 016802. 
[76] Prus O., Yaish Y., Reznikov M., Sivan U., and Pudalov V., Phys. Rev. B, 67 

(2003) 205407. 

[77] Lakner M., v.Lohneysen H., Langenfeld A., and Wolfle P., Phys. Rev. B, 50 
(1994) 17064. 

[78] Castellani C, Di Castro C, Forgacs C, and Sorella S., Solid State 

Communications, 52 (1984) 261. 
[79] Kopietz P., Phys. Rev. Lett, 81 (1998) 2120. 



